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Summary. — In the last two decades the search for neutrinoless double beta decay 
has evolved into one of the highest priorities for understanding neutrinos and the 
origin of mass. The main reason for this paradigm shift has been the discovery of 
neutrino oscillations, which clearly established the existence of massive neutrinos. 
An additional motivation for conducting such searches comes from the existence of 
an unconfirmed, but not refuted, claim of evidence for neutrinoless double decay in 
^^Ge. As a consequence, a new generation of experiments, employing different detec- 
tion techniques and /3/3 isotopes, is being actively promoted by experimental groups 
across the world. In addition, nuclear theorists are making remarkable progress in 
the calculation of the neutrinoless double beta decay nuclear matrix elements, thus 
eliminating a substantial part of the theoretical uncertainties affecting the particle 
physics interpretation of this process. In this report, we review the main aspects 
of the double beta decay process and some of the most relevant experiments. The 
picture that emerges is one where searching for neutrinoless double beta decay is rec- 
ognized to have both far-reaching theoretical implications and promising prospects 
for experimental observation in the near future. 

PACS 23.40. -s - /3 decay; double /3 decay; electron and muon capture. 
PACS 14 . 60 . Pq ~ Neutrino mass and mixing. 
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1. Introduction 



Neutrinoless double beta decay (/3/30z^) is a hypothetical nuclear transition in which 
two neutrons undergo /3-decay simultaneously and without the emission of neutrinos. If 
realized in Nature, this transition would be extremely rare: the most constraining lower 
bound on the neutrinoless double beta decay half-life, in ^^Gc, is T"^2 > 1-9 x 10^^ years 
(at 90% confidence level, from [1]). The importance of /3/30i^ searches goes beyond its 
intrinsic interest, as it is the only practical way to reveal experimentally that neutrinos 
are Majorana particles. If i/ is a field describing a neutrino, stating that the neutrino 
is a Majorana particle is equivalent to saying that the charge-conjugated field — that 
is, a field with all charges reversed — also describes the same particle: v = ly'^. If such 
Majorana condition is not fulfilled, we speak instead of Dirac neutrinos. 
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The theoretical impUcations of experimentaUy estabhshing (3(30iy would be profound. 
In a broad sense, Majorana neutrinos would constitute a new form of matter, given 
that no Majorana fermions have been observed so far. Also, /3/30;/ observation would 
prove that total lepton number is not conserved in physical phenomena, a fact that could 
be linked to the cosmic asymmetry between matter and antimatter. Finally, Majorana 
neutrinos would mean that a new physics scale must exist and is accessible in an indirect 
way through neutrino masses. 

In addition to theoretical prejudice in favor of Majorana neutrinos, there are other 
reasons to hope that experimental observation of Pj3Qv is at hand. Neutrinos are now 
known to be massive particles, thanks to neutrino oscillation experiments. If PfSOv 
is mediated by the standard light Majorana neutrino exchange mechanism, a non-zero 
neutrino mass would almost certainly translate into a non-zero /3/30i^ rate. While neutrino 
oscillation phenomenology is not enough per se to provide a firm prediction for what such 
a rate should be, it does give us hope that a sufficiently fast one to be observable may 
be realized in Nature. Furthermore, /3/30i^ may have been observed already: there is 
an extremely intriguing, albeit controversial, claim for PfSQv observation in ^^Ge that is 
awaiting unambiguous confirmation by future PPOv experiments. 

The profound theoretical implications of massive Majorana neutrinos, and the pos- 
sibility that an experimental observation is at hand, has triggered a new generation of 
jSpOv experiments. At the time of writing this report, this new generation of experi- 
ments spans at least half a dozen isotopes, and an equally rich selection of experimental 
techniques, ranging from the well-established germanium calorimeters, to xenon time 
projection chambers. Some of the experiments are already running or will run very soon. 
Some of them are still in their R&D period. Some of them push to the limit the technique 
they use, in particular concerning the target mass. Others are easier to scale up. All of 
them claim to be sensitive to very light neutrino masses, by assuming that they can do 
one to three orders of magnitude better in background suppression and by significantly 
increasing their target mass, compared to previous experiments. In this report we review 
the state-of-the-art of this exciting and rapidly changing field. 

This review is organized as follows. The introductory material is covered in sects. 2 
and 3. The key particle physics concepts involving massive Majorana neutrinos and 
neutrinoless double beta decay are laid out here. The current experimental knowledge 
on neutrino masses, lepton number violating processes in general, and jSpOv in particular, 
is also described in sects. 2 and 3. Sections 4, 5 and 6 cover more advanced topics. The 
theoretical aspects of the nuclear physics of j3j3Qv are discussed in sect. 4. Sections 5 
and 6 deal with experimental aspects of fS/SOv, and can be read without knowledge of 
sect. 4. An attempt at a pedagogical discussion of experimental ingredients affecting 
PfSOiy searches is made in sect. 5. Section 6 adds a description of selected new- generation 
experimental proposals, together with a comparison of their physics case. 

2. Massive neutrinos 

2'1. Current knowledge of neutrino mass and mixing. - Neutrinos are the lightest 
known elementary fermions. Neutrinos do not carry any electrical charge, do not undergo 
strong interactions, and are observable only via weak interactions. In the Standard Model 
of elementary particles, neutrinos are paired with charged leptons in weak isodoublets. 
Experimentally, we know that only three light (that is, of mass < mz/2, where mz is 
the Z boson mass) active neutrino families exist. 

More recently, neutrino oscillation experiments have unambiguously demonstrated 
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that neutrinos are massive particles (see, for example, [2]). Because of the interferometric 
nature of neutrino oscillations, such experiments can only measure neutrino squared mass 
differences and not the absolute neutrino mass scales. Solar and reactor experiments 
have measured one mass splitting, the so-called solar mass splitting, to be: Am^^^ = 
1712 ~ "^1 — i'^-'^^-o 'ie) ^ 10^^ eV^. Atmospheric and accelerator-based experiments 
have measured a different mass splitting, the so-called atmospheric mass splitting, to be: 
lAm^^^l = |m| - {ml + ml)/2\ = (2.35lj] J^) x fO'^ eV^ > Am^^j. In the standard 
3-neutrino oscillations paradigm, those are the only two independent mass splittings 
available. The best-fit values and la ranges quoted were obtained from a recent global 
3-neutrino fit [3]. 

The observation of neutrino flavor oscillations also imply that the neutrino states 
participating in the weak interactions {flavor eigenstates) are different from the neutrino 
states controlling free particle evolution {mass eigenstates). In other words, the three 
weak eigenstates \va), a — e, ^, r, can be expressed as a linear combination of the three 
mass eigenstates i = 1,2,3: 



where ?7 is a 3 x 3, unitary, neutrino mixing matrix, that is different from unity. Equation 
1 implies the violation of the individual lepton flavors La , but not necessarily the violation 
of total lepton number L = La = Lf,+ L^^ -\- L^. The 3x3 neutrino mixing matrix is 
usually parametrized in terms of 3 Euler angles (1^12, i^is, 1^23) and 3 phases {6, a2i, a^i) 
(see, for example, [4]). If the massive neutrinos are Dirac particles (see sect. 2'2), only 
the Dirac phase 6 is physical and can be responsible for CP violation in the lepton sector. 
If the massive neutrinos are Majorana particles (sect. 2'2), the two additional Majorana 
phases (q!2i, aai) arc also potentially observable. 

Neutrino oscillation experiments have measured with reasonably good accuracy the 
flavor content of the neutrino mass states participating in 3-ncutrino mixing. Atmo- 
spheric and accelerator-based neutrino oscillation experiments are mostly consistent with 

— )■ v-r oscillations only (see, however, [5]). They have therefore measured the muon 
flavor content of the 1^3 mass state to be |J7^3p — 0.5, and that such mass state has 
little (if non-zero) electron flavor content, |J7e3p — 0. On the other hand, solar neutrino 
oscillation experiments are consistent with i^e — >■ and/or ~> J^t oscillations. They 
have measured |J7e2p — 1/3. The remaining elements of the leptonic mixing matrix can 
approximately be derived, given (|C/p3|, \Ue3\, |C^e2|), assuming unitarity. More precisely, 
according to [3], the best-fit values and la ranges in the neutrino mixing parameters mea- 
sured via neutrino oscillations are: pesl'^ = 0.025 ±0.07, |C/^3|V(1 - Pesl'^) = 0.42^°;°^ 
and |?7e2|V(l " |C^e3p) = 0.312l°:oi[6- The value of the Dirac CP-violating phase 6, also 
potentially observable in neutrino oscillation experiments, is currently unknown. 

The current knowledge on neutrino masses and mixings provided by neutrino oscil- 
lation experiments is summarized in fig. 1. The diagram shows the two possible mass 
orderings that arc compatible with neutrino oscillation data, with increasing neutrino 
masses from bottom to top. In addition, the electron, muon and tau flavor content of 
each mass eigenstate is also shown, according to the best-fit values in reference [3]. 

To complete our knowledge on neutrino masses, two pieces of information remain to 
be known: the neutrino mass ordering and the absolute value of the lightest neutrino 
mass. 



(1) 
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Fig. 1. - Knowledge on neutrino masses and mixings from neutrino oscillation experiments. 
Panels (a) and (b) show the normal and inverted mass orderings, respectively. Neutrino masses 
increase from bottom to top. The electron, muon and tau flavor content of each neutrino mass 
eigenstate is shown via the red, green and blue fractions, respectively. 



Concerning the neutrino mass ordering, current neutrino oscillation results cannot 
differentiate between two possibilities, usually referred to as normal and inverted or- 
derings (see fig. 1). In the former, the gap between the two lightest mass eigenstates 
corresponds to the small mass difference, indicated by solar experiments (Am^Qj), while 
in the second case the gap between the two lightest states corresponds to the large mass 
difference, inferred from atmospheric experiments (Am^^jjj). While we do not know at 
present whether vt, is heavier or lighter than vi, we do know that V2 is heavier than 
vi, thanks to matter effects affecting the propagation of neutrinos inside the Sun. The 
exploitation of the same type of matter effect on future accelerator-based neutrino exper- 
iments may allow us to experimentally establish the neutrino mass ordering in the future, 
provided that |J7e3| > 0. In the particular case in which the neutrino mass differences 
are very small compared with its absolute scale, we speak of the degenerate spectrum. 

The absolute value of the neutrino mass scale can instead be probed via neutrinolcss 
double beta decay searches, cosmological observations and beta decay experiments. Only 
upper bounds on the neutrino mass, of order '--^l eV, currently exist. Constraints on the 
lightest neutrino mass coming from neutrinoless double beta decay will be discussed in 
sect. 3'3. In the following, we briefly summarize cosmological and beta decay constraints. 

Primordial neutrinos have a profound impact on cosmology since they affect both 
the expansion history of the Universe and the growth of perturbations (see, for instance, 
reference [6]). Cosmological observations can probe the sum of the three neutrino masses: 



(2) ^cosmo — ^ ^ 



mi 
1 



Cosmological data are currently compatible with massless neutrinos. Several upper 
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(a) (b) 

Fig. 2. - Constraints on the lightest neutrino mass miight coming from a) cosmological and b) jS 
decay experiments. The red and green bands correspond to the normal and inverted orderings, 
respectively. The rricosmo upper bound in panel (a) is from [7] , and translates into a mught upper 
limit shown via the vertical band in the same panel. The cosmological constraint on mught is 
also shown in panel (b), together with the upper limit on mg from tritium /? decay experiments 
[4]. 



limit values on mcosmo can be found in the literature, depending on the details of the 
cosmological datasets and of the cosmological model that were used in the analysis. A 
conservative upper limit on mcosmo of 1.3 cV at 95% confidence level [7] is obtained 
when CMB measurements from the Wilkinson Microwave Anisotropy Probe (WMAP) 
are combined with measurements of the distribution of galaxies (SDSSII-BAO) and of 
the Hubble constant iJo (HST), in the framework of a cold dark matter model with dark 
energy whose equation of state is allowed to differ from —1. The relationship between 
TOcosmo, defined in eq. (2), and the lightest neutrino mass mught — that is, mi (ms) in 
the case of normal (inverted) ordering — is shown in fig. 2a. The two bands correspond 
to the normal and inverted orderings, respectively. The width of the bands is given by 
the 3(7 ranges in the mass oscillation parameters Am^^j and Am^^j^^ [3]. The horizontal 
band in fig. 2(a) is the upper limit on mcosmo- In this quasi-degenerate regime, this upper 
bound implies that mught — mcosmo/3 ^ 0.43 eV at 95% CL, as shown by the vertical 
band in fig. 2(a). 

The neutrino mass scale can also be probed in laboratory-based experiments (see, for 
example, [8]). The differential electron energy spectrum in nuclear j3 decay experiments 
is affected both by the neutrino masses and by the mixings defining the electron neutrino 
state in terms of mass eigenstates. In this case, the mass combination probed is given 
by: 

3 

(3) m|^^|C/,.|V 

1=1 

The relationship between in eq. (3) and m^ght is shown in fig. 2(b). Again, the 
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results of a recent global fit to neutrino oscillation data [3] are used to determine the Scr 
bands for both the normal and inverted orderings. From the experimental point of view, 
the region of interest for the study of neutrino properties is located near the /? endpoint. 
The most sensitive searches conducted so far are based upon the decay of tritium, via 
— He~''e~Pe, niostly because of the very low /3 endpoint energy of this element (18.6 
keV). As for cosmology, /3 decay searches of neutrino mass have so far yielded negative 
results. The horizontal band in fig. 2(b) comes from the two current best limits from 
the Troitsk [9] and Mainz [10] experiments. The combined limit is < 2 eV at 95% 
CL [4]. The resulting constraint on miight is less stringent than the cosmological one. 
The KATRIN experiment [8] should be able to improve the (and therefore mught) 
sensitivity by roughly an order of magnitude in the forthcoming years, thanks to its 
better statistics, energy resolution, and background rejection. 

2'2. The origin of neutrino mass: Dirac versus Majorana neutrinos . ~ Are neutrinos 
their own antiparticles? If the answer is no, we speak of Dirac neutrinos. If the answer 
is yes, we speak of Majorana neutrinos. Both possibilities exist for the neutrino, being 
electrically neutral and not carrying any other charge-like conserved quantum number. 
Whether neutrinos are Majorana or Dirac particles depends on the nature of the physics 
that give them mass, given that the two characters are physically indistinguishable for 
massless neutrinos. 

In the Standard Model, only the negative chirality component ^'l of a fermion field 
vf" = 'iiji + is involved in the weak interactions. A negative (positive) chirality field 
*L(fi) is a field that obeys the relations Pl{r)'^l{r) = *L(it) and Pr{l)'^ l(r) = 0, where 
Pl ~ {\ — 75)/2 and Pn = (1 + 75)/2 arc the positive and negative chiral projection 
operators. 

For massless neutrinos (see, for example, [11]), only the negative chirality neutrino 
field vl is needed in the theory, regardless of the Dirac/Majorana nature of the neutrino 
discussed below, since neutrinos only participate in the weak interactions. This field 
describes negative helicity neutrino states {vl) and positive helicity antincutrino states(^). 
The positive and negative helicity states are eigenstates of the helicity operator h = a ■ p 
with eigenvalues ±1/2, respectively, where a is the neutrino spin and p the neutrino 
momentum direction. The fact that vl annihilates particles of negative helicity, and 
creates antiparticles with positive helicity, is not inconsistent with Lorentz invariance, 
given that the helicity is the same in any reference frame for a fermion travelling at the 
speed of light. In the Standard Model with massless neutrinos, positive helicity neutrinos 
and negative helicity antineutrinos do not exist. As a consequence, and since a negative 
helicity state transforms into a positive helicity state under the parity transformation, 
the chiral nature of the weak interaction (differentiating negative from positive chirality) 
implies that parity is maximally violated in the weak interactions. 

For relativistic neutrinos of non-zero mass m, the neutrino field participating in the 
weak interactions has still negative chirality, vl , but there are sub- leading corrections to 
the particle annihilation/creation rules described above. The state {vl) that is annihi- 
lated by the negative chirality field i'l is now a linear superposition of the —1/2 and -t-1/2 



(^) As customarily done, we use the subscript "L" to denote both negative helicity states 
[z^l) and negative cliirality fields vl, since the terms left-handed helicity states and left-handed 
chirality fields are also commonly used. Similarly, we denote positive helicity states and positive 
chirality fields with the subscript "R" , as in "right-handed" . 
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helicity states. The +1/2 helicity state enters into the superposition with a coefficient 
cx m/ where E is the neutrino energy, and is therefore highly suppressed. 

Neutrino mass terms can be added to the Standard Model Lagrangian in two ways 
(see, for example, [12]). The first way is in direct analogy to the Dirac masses of quarks 
and charged leptons, by adding the positive chirality component of the Dirac neu- 
trino field, describing predominantly positive helicity neutrino states and predominantly 
negative helicity antineutrino states that do not participate in the weak interactions: 

(4) -Cd =mD {i^VR + TTJ^VL ) , 

where mu = yv/\/2 > 0, j/ is a dimensionless Yukawa coupling coefficient and v/V^ is the 
vacuum expectation value of the neutral Higgs field after electroweak symmetry breaking. 
In eq. (4), and vr are, respectively, the negative and positive chirality components of 
the neutrino field v. The chiral spinors vl and iTi have only two independent components 
each, leading to the four independent components in the spinor ly. This is different from 
the the case of massless neutrinos, where only the 2-component spinor i^l was needed. 

The second way in which neutrino mass terms can be added to the Standard Model 
Lagrangian is unique to neutrinos. Majorana first realized [13] that, for neutral particles, 
one can remove two of the four degrees of freedom in a massive Dirac spinor by imposing 
the Majorana condition: 



where i^''^ ~ Cv^ ~ C{^'^)^v* is the CP conjugate of the field C is the charge- 
conjugation operator, and {vlY H'^rT) ^i^s positive (negative) chirality. This result can 
be obtained by decomposing both the left-hand and right-hand sides of eq. (5) into their 
chiral components, yielding: 



and therefore proving that the positive chirality component of the Majorana neutrino 
field is not independent of, but obtained from, its negative chirality counterpart ly^. 
By substituting eq. (6) into the mass term in eq. (4), we obtain a Majorana mass term: 



where ttil is a free parameter with dimensions of mass. This lagrangian mass term 
implies the existence of a weak isospin triplet scalar (a Higgs triplet), with a neutral 
component acquiring a non-vanishing vaccum expectation value after electroweak sym- 
metry breaking. Equation (7) represents a mass term constructed from negative chirality 
neutrino fields alone, and we therefore call it a negative chirality Majorana mass term. 
If positive chirality fields also exist and are independent of negative chirality ones, this 
is not the only possibility. In this case, we may also construct a second Majorana mass 
term, a positive chirality Majorana mass term: 



(5) 



(6) 



(7) 




(8) 
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Fig. 3. - The difference between Dirac (left) and Majorana (right) massive neutrinos in a scat- 
tering experiment. See text for details. Adapted from [14]. 



In the Standard Model, right-handed fcrmion fields such as vr are weak isosopin 
singlets. As a consequence, and in contrast with m^) or m/,, the mass parameter is 
therefore not connected to a Higgs vaccum expectation value, and could be arbitrarily 
high. All three mass term in eqs. (4), (7) and (8) convert negative chirality states into 
positive chirality ones(^). Chirality is therefore not a conserved quantity, regardless of the 
Dirac/Majorana nature of neutrinos. Furthermore, the Majorana mass terms in eqs. (7) 
and (8) convert particles into their own antiparticles. As stated previously, they are 
therefore forbidden for all electrically charged fermions because of charge conservation. 
But not only: processes involving Majorana mass terms violate the Standard Model total 
lepton number L = Le + -^/^ + -^t by two units (| AL| = 2), which is not a good quantum 
number anymore. 

Which of the mass terms allowed in theory, among Co-, Cl and Cr in eqs. (4), (7) 
and (8) exist in nature? What are the numerical values of the corresponding coupling 
constants m/^, m^,, mi?? These questions can in principle be answered experimentally. 
Majorana and Dirac massive neutrinos will in fact have different Standard Model inter- 
actions. Let us consider for now an instructive, albeit unrealistic, scattering experiment, 
see fig. 3. 

In the Dirac case. Standard Model interactions conserve lepton number L, with 
L{v) = L{1~) = —L{v) = —L{1^) = -fl, where indicates charged leptons. Particles 
are then identified as neutrinos or antineutrinos in accordance with the process through 
which they are produced. Charged-current interactions of Dirac neutrinos (as opposed 
to antineutrinos) produce only and carry a well-defined lepton number L = — 1, 
and viceversa. As shown in fig. 3, for Dirac neutrinos we would thus have four mass- 
degenerate states: for each of the two available helicity states(^), two distinct parti- 
cle/antiparticle states characterized by a different L value would be available. Standard 
Model interactions of neutrino (as opposed to antineutrino) states of positive helicity 
would have, however, much weaker /"-producing interactions with matter compared to 
neutrino states of negative helicity, as indicated by the coefficients in fig. 3. On the other 
hand, we have seen that in the Majorana case L is not conserved. We would only have 
two mass-degenerate states, defined by the two available helicity states, see fig. 3. 



(^) This is because the charge conjugate of a field with a given chirality, such as the ones 
appearing in eqs. 7 and 8, always has the opposite chirality. 

(^) As mentioned above, the weak interaction is maximally parity violating, therefore the two 
helicity states are distinguishable. 
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Fig. 4. - Hierarchical structure of fermion masses. Only upper bounds for neutrino masses exist. 
The figure assumes a normal ordering for neutrino masses. Values taken from [4] . 



Given these difFcrences between Dirac and Majorana massive neutrinos, can we es- 
tablish which of the two possibilities is realized in Nature via a scattering experiment? 
In practice, no. The reason is that l~ production from positive helieity Dirae neutrinos 
(and l'^ production from negative helieity Dirac antineutrinos) is expected to be highly 
suppressed in the ultra-relativisitc limit, and cannot be experimentally observed. Ex- 
perimentally, all we know is that the neutral particle produced in association with a 
l'^ produces, when interacting, a l~ . In the Dirac case, lepton number conservation is 
assumed and such neutral particle is identified as the neutrino, with L = —I. In the 
Majorana case, such neutral particle is instead identified as the negative helieity state, 
interacting differently from its positive helieity counterpart. Both interpretations are 
viable, and what happens when a neutrino interacts can be understood without invoking 
a conserved lepton number [15]. 

2'3. The see-saw mechanism. - Neutrino masses, although not measured yet, are 
known to be small, of the order of 1 eV or less, see sect. 2'1. Such mass values are much 
smaller than the masses of all other elementary fermions, see fig. 4. The explanation of 
neutrino masses via Dirac mass terms alone require neutrino Yukawa couplings of the or- 
der of 10~^^ or less. The current theoretical prejudice is that neutrino Yukawa couplings 
with <C 1 and <^ yi arc unnatural, if not unlikely. 

The so-called see-saw mechanism provides a way to accommodate neutrino masses 
that is considered more natural. The simplest realization of the see-saw model is to add 
both a Dirac mass term and a positive chirality mass term to the Lagrangian, as given 
by eqs. (4) and (8), respectively, for each of the three neutrino flavors. This is sometimes 
called the type I see-saw mechanism, where we take rriL = 0, m^) 7^ 0, and mu 7^ 0. In 
this case, the neutrino mass terms can be recast in the matrix form: 



(9) 



-Cd+k = 2 (^l)' M Ml + h.c. , 
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where the matrix M has the form: 

(10) M = 
and the negative chirahty vector Ml is 

(11) Ml = 



Tod 
Tod TOfl 



VL 



The chiral fields vl and i/^ do not have a definite mass, since they are coupled by 
the Dirac mass term. In order to find the fields vil and iVi^ with definite masses toi 
and Ml, respectively, it is necessary to diagonalize the mass matrix in eq. (9). In other 
words, it is necessary to find a unitary mixing matrix U such that: 

(12) (7 
where: 

(13) J\fL=l^nL, with ul 



^IL 

Nil 



For each neutrino flavor, two fields of definite chirality and definite mass are therefore 
obtained, and the diagonalized mass terms can be written as: 

1 



(14) -Cd+r = 2 [mi {i^ilY- >^il + Mi {NilY Nil) + h.c. , 

Both terms in eq. (14) have the same form as the pure negative chirality Majorana 
mass term in eq. (7). In other words, both mass eigenfields vil and iVi^ are equal to 
their CP-conjugate fields, and thus both describe Majorana particles. The insertion of 
a Dirac mass term and a positive chirality Majorana mass term in the Lagrangian for 
massive neutrinos has resulted in Majorana particles. 

Since the positive chirality fields are electroweak singlets in the Standard Model, the 
Majorana mass of the neutrino described by such field, to^j, may be orders of magnitude 
larger than the electroweak scale. In the so-called see-saw limit, we assume that neutrino 
Yukawa couplings are of the order of the charged fcrmion couplings, and that toj^ ^ niD 
is of the order of some high mass scale where new physics responsible for neutrino masses 
is supposed to reside. In this approximation, the see-saw mechanism naturally yields 
a small mass eigenvalue mi ~ mj-,/mii <^ niD for a predominantly negative helicity 
neutrino mass state, and a large mass eigenvalue Mi mfj for a predominantly positive 
helicity (and therefore sterile) neutrino mass state. A very heavy A'^i corresponds to a 
very light i>i and viceversa, as in a see-saw. 

The see-saw mechanism presented above can easily be generalized from the one-family 
case that we discussed to three neutrino species, yielding the three light neutrinos Vi we 
are familiar with, and three heavy neutrinos Ni, with i = 1,2,3 [12]. In this case, the 
neutrino mass matrix in eq. 9 is a 6 x 6 mass matrix of the form: 



(15) M 



(A/D)^ 
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where and are now 3x3 complex matrices, and the 6-component vector of 
negative chirality fields has the form: 



(16) Nl = [ ^^"^y ) , with VL =\yt.L\ and (vrY = | {v.^rY 




{^s.rY 

S'^ssrT 



In eq. 16, the subscripts e, /i, r label the active neutrino flavors, while the subscripts 
■Si, 32, S3 indicate sterile states that do not participate in the weak interactions. The 
mass matrix is diagonalized via a 6 x 6 mixing matrix V analogous to U in eq. 12, 
where the three negative chirality fields and the three positive chirality fields arc now 
expressed in terms of the negative chirality components of 6 massive neutrino fields 
i^iL, I = 1, . . . , 6. In the see-saw limit where the eigenvalues of are much larger than 
those of M^, the 6x6 mass matrix in eq. 15 can be written in block-diagonal form 
M ~ diag(Miight, Afhoavy), where the two 3x3 mass matrices of the light and heavy 
neutrino sectors are practically decoupled, and given by Mught — —{AI^)'^{M^)~^M^ 
and Mhoavy — M^, respectively. 

For the low-energy phenomenology, it is sufficient to consider only Mught, sometimes 
called the neutrino mass matrix m^, that is the 3x3 matrix in the flavor basis which is 
diagonalized by the matrix U : 

(17) Afiight U = diag(TOi, 7712, ma) , 

where the neutrino mixing matrix U appearing in eq. (17) is the same matrix defined in 
eq. (1), and toi, m2, are three light neutrino mass eigenvalues discussed in sect. 2'1. 

An important assumption in the simplest realization of the see-saw mechanism de- 
scribed above is that = 0. This assumption is not arbitrary, and directly follows 
from enforcing the gauge symmetries of the Standard Model, see for example [12]. In 
models with a left-right symmetric particle content, this type I sce-saw mechanism is 
often generalized to a type II see-saw, where an additional direct mass term tul for the 
light neutrinos is present. 

2'4. Leptogenesis . - Inflationary models of the Universe predict matter and antimatter 
to be equally abundant at the very hot beginning, given that any potential initial asym- 
metry would have been diluted away by inflation. However, the observable Universe today 
is almost entirely made of matter! This matter dominance today is consistent with the 
small level of baryon asymmetry that is inferred from BBN and CMB observations, given 
that annihilation of matter with antimatter would have left us in a matter-dominated 
Universe today. The baryon asymmetry has been precisely measured: 

(18) r, = ^ 273.9 • IQ-^^flbh^ ^ (6.19 ± 0.15)10-1° 

where ng, ng, arc the number densities of baryons, antibaryons and photons, flf, ~ 
(0.0458 ± 0.0016) is the fraction of the critical energy density carried by baryons, and 
h = iJo/100 km • s-i • Mpc'^ = (0.742 ± 0.036) is the Hubble parameter, where Hq is 
the Hubble constant today. 

What caused this matter-antimatter asymmetry in the early Universe? The baryon 
asymmetry could have been induced by a lepton asymmetry: leptogenesis (see, for ex- 
ample, [16, 17]). If neutrinos are Majorana particles, the decays of the heavy Majorana 
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(a) (b) (c) 

Fig. 5. - Feynman diagrams that contribute to the lepton number asymmetry through the 
decays of the heavy Majorana neutrino A^i into the Higgs (j) plus leptons la- The asymmetry is 
generated via the interference of the tree-level diagram (a) with the one-loop vertex correction 
(b) and the self-energy (c) diagrams (adapted from [16]). 



neutrinos into leptons la plus Higgs particles (p in the early Universe provides an ideal 
scenario for leptogenesis. Heavy Majorana neutrinos are their own anti-particles, so they 
can decay to both and la4' final states. If there is an asymmetry in the two decay 
rates, a net lepton asymmetry will be produced. Figure 5 shows the processes that would 
contribute to a net lepton asymmetry in the presence of heavy Majorana neutrino de- 
cays, in the simplest case where the asymmetry is dominated by the decay of the lightest 
among the three heavy neutrinos, A'^i. Finally, this lepton asymmetry can be efficiently 
converted into a baryon asymmetry via the so-called sphaleron processes (see [16, 17] for 
details). 

In more detail, for leptogenesis to occur, three conditions must be met. These con- 
ditions directly follow from the ingredients that are required to dynamically generate a 
baryon asymmetry {Sakharov's conditions [18]): 

1. Presence of lepton number violating processes; 

2. Beyond-SM sources of CP violation ("*); 

3. Departure from thermal equilibrium. 

The decay of heavy Majorana neutrinos can provide all of these conditions, namely 

1. Total lepton number is violated in these decays; 

2. CP can be violated in these decays, provided that there is more than one heavy 
Majorana field; 

3. Departure from thermal equilibrium is obtained if the decay rate is slower than the 
expansion rate of the Universe at the time of decoupling, occurring for T Mi, 
where T is the temperature of the Universe's thermal bath, and Mi is the mass of 
the lightest among the three heavy neutrinos. 

In order to be fully successful, any theory of leptogenesis must be able to explain the 
observed magnitude of baryon asymmetry given in eq. (18). Leptogenesis via heavy 



(*) CP violation is allowed by the Standard Model and has been measured; however, the magni- 
tude of such CP-violating effects is far too small to provide the necessary amount of leptogenesis. 
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Majorana neutrino decays is in principle able to do this. In this case, the asymmetry in 
lepton flavor a produced in the decay of A'^i, defined as: 

_ r(jVi ^ (bL) - r(iVi ^ 4>L) 

should be of order \£aa\ > 10^^ [17], where the factors F in eq. (19) stand for the decay 
rates into the corresponding iVi decay final states. It is at present unclear whether there 
is a direct connection between the high-energy CP-violating processes responsible for the 
asymmetry in the early Universe of eq. (19), and the low-energy CP-violating processes 
that may potentially affect laboratory-based experiments. Nonetheless, the discovery 
of CP violation in the lepton sector via neutrino oscillations on the one hand, and the 
discovery of the Majorana nature of neutrinos via ncutrinoless double beta decay on the 
other, would undoubtedly strengthen the case for leptogenesis as a source of the baryon 
asymmetry of the Universe. 

2'5. Lepton number violating processes. - We have seen that Majorana mass terms 
induce lepton number violating processes of the type |AL| = 2. The heavy neutrino 
decay needed for leptogenesis, discussed in sect. 2 '4, requires Majorana mass terms and 
is therefore an example of a lepton number violating process. However, heavy neutrino 
decay is unobscrvable in a laboratory-based experiment, given the tremendous energies 
needed for heavy neutrino production. A number of more promising lepton number 
violating processes have been proposed to probe the Majorana nature of neutrinos. The 
best known example is fi^Qv, subject of this review and introduced in sect. 3. We 
anticipate that PPOv is considered the most promising probe of the Majorana nature of 
neutrinos. However, and because of neutrino mixing, the phenomenology associated with 
|Ai| = 2 processes is very rich. The basic process with |AL| = 2 is mediated by [19]: 

(20) W-W~ ^ l-lj 

and we can categorize such processes according to the lepton flavors (a, /3) involved. 
Assuming no lepton flavor violating contributions other than light Majorana neutrino 
exchange, the matrix element for the generic |AL| = 2 process in eq. (20) is proportional 
to the element (a,/3) of the neutrino mass matrix: 

3 

(21) (m.)„^ ^ ([/*diag(mi,m2,m3)C/t)^^ ^Y.^:^U;,m, 

1=1 

where ruu = M''^''* is the matrix appearing in eq. (17), Uai are the elements of the 3x3 
neutrino mixing matrix appearing in eq. 1, and are the three light neutrino masses. 
In a sense, this effective neutrino mass definition provides a metric to compare the sensi- 
tivity of various |AL| = 2 processes. The processes with the most competive constraints 
on |Ai| = 2 processes involving the flavors (a,/3) are reported in table I. 

As is apparent in table I, indeed the constraint on the effective Majorana mass rriee 
coming from /3P0iy searches outperforms by several orders of magnitude other searches 
involving a different flavor combination (a,/3). The most important reason behind this 
is of statistical nature. While it is possible to amass macroscopic quantities of a /3f3 
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Table I. - Current bounds on effective neutrino masses from total lepton number violating 
processes, organized according to the flavors involved. Numbers taken (or derived) from [4, 19]. 



Flavors 


Exp. technique 


Exp. bound 


Mass bound (eV) 


(e,e) 




Ti/zC^'Ge ^ '"'Se + 26") > 1.9 x 10^^ yr 


\mee \ < 3.6 X 10"^ 


(e,M) 


— >■ e'^ conversion 


r(Ti + fi- ^e+ + Gag,) / 
r(Ti + ^"capture) < 1.7 x 10"^^ 


|m,^| < 1.7 X 10^ 


(e,r) 


Rare r decays 


r{T- e+7r-7r~)/rtot < 8.8 x 10"* 


|mer| < 2.6 X 10^^ 




Rare kaon decays 


r{K+ ^ 7r-M+M+)/rtot < 1.1 X 10-9 


|m^^| < 2.9 X 10* 




Rare t decays 


r(r- ^ /i+7r-7r-)/rtot < 3.7 x 10"* 


|mer| < 2.1 X 10^^ 


{r,r) 


none 


none 


none 



emitter to study (SpOiy decay (as we will see, even a ton of isotope is in the cards of 
several experiments), this is not the case for the other experimental techniques listed in 
table I. Nevertheless, it is important to keep exploring lepton flavor violating processes 
other than fS/SOi/ for two reasons. First, it is in principle possible that phase cancellations 
are such that <^ rriap with (a,/3) 7^ (e, e), making the search for /3/30h' much less 
favorable than others because of Nature's choice of neutrino masses and mixings. Second, 
this effort may possibly lead to the identification of an even most promising experimental 
probe of lepton flavor violation in the future. 

3. Neutrinoless double beta decay 

3'1. Double beta decay modes. - Double beta decay is a rare nuclear transition in which 
a nucleus with Z protons decays into a nucleus with Z + 2 protons and the same mass 
number A. The decay can occur only if the initial nucleus is less bound than the final 
nucleus, and both more than the intermediate one, as shown in fig. 6. Such a condition 
is fulfilled by 35 nuclides in nature because of the nuclear pairing force (see sect. 4'2.2), 
ensuring that nuclei with even Z and N are more bound than the odd-odd nuclei with 
the same A = Z + N. 

The standard decay mode (/3/32i/), consisting in two simultaneous beta decays, 

(22) {Z,A)^ {Z + 2,A) + 2 e- + 2Ve, 

was first considered by Maria Goeppert-Mayer in 1935 [20]. Total lepton number is 
conserved in this mode, and the process is allowed in the Standard Model of particle 
physics. This process was first detected in 1950 using gcochemical techniques [21]. The 
first direct observation of j3p2v events, in ®^Se and using a time projection chamber, 
did not happen until 1987 [22]. Since then, it has been repeatedly observed in several 
nuchdes, such as "^^Ge, ^°°Mo or ^^°Nd. Typical hfetimes are of the order of 10^^-10^° 
years, the longest ever observed among radioactive decay processes. For a list of I3(i2v 
half-lives measured in several isotopes, see table II [23]. The longest half-life in tab. II is 
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Fig. 6. - Atomic masses of A — 136 isotopes. Masses are given as differences witfi respect to the 
most bound isotope, ^'^^Ba. The red (green) levels indicate odd-odd (even-even) nuclei. The 
arrows P~ P~ indicate nuclear decays accompanied by electron, positron and double 

electron emission, respectively. The arrows EC indicate electron capture transitions. 



Table II. - Current best direct measurements of the half-Ufe of I3j32v processes. The values 
reported are taken from the averaging procedure described m [23j. 



Isotope 


T^h (year) 


Experiments 


*»Ca 
^«Ge 
8=Se 
36 Zr 
iooMo 
i"Cd 

130 

"6Xe 


(4.4t0;6) X 10^3 
(1.5 ±0.1) X 10^^ 
(0.92 ± 0.07) X 10^° 
(2.3 ±0.2) X 10^3 
(7.1 ±0.4) X 10^* 
(2.8 ±0.2) X 10^3 
(6.8ti:?) X 10^° 
(2.11 ±0.21) X 10^^ 
(8.2 ±0.9) X 10^^ 


Irvine TPC [28], TGV [29], NEM03 [30] 
PNL-USC-ITEP-YPI [31], IGEX [32], H-M [33] 
NEM03 [34], Irvine TPC [35], NEM02 [36] 

NEM02 [37], NEM03 [38] 
NEM03 [34], NEMO-2 [39], Irvine TPC [40] 
NEM03 [30], ELEGANT [41], Solotvina [42], NEM02 [43] 
CUORICINO [44], NEM03 [45] 

EXO-200 [24] 
Irvine TPC [40], NEM03 [46] 



the one for ^^^Xe, whicli has been measured for the first time only in 2011 [24] (^). 
The neutrinoless mode (/3/30z^), 

(23) (Z,A) ^ (Z±2,A)±2 e", 

was first proposed by W. H. Furry in 1939 [47] as a method to test Majorana's theory 



(^) The 10% accuracy in the ^'^^Xe /3/32i/ decay measured half-life in [24] should be contrasted 
with a spread of more than one order of magnitude in the corresponding theoretical expectations 
from several nuclear structure calculations [25, 26, 27]. 
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[13] applied to neutrinos. In contrast to the two-neutrino mode, the neutrinoless mode 
violates total lepton number conservation and is therefore forbidden in the Standard 
Model. Its existence is linked to that of Ma jorana neutrinos (see sect. 3'2). No convincing 
experimental evidence of the decay exists to date (see sect. 3'5). 

The two modes of the /3/3 decay have some common and some distinct features [48]. 
The common features are: 

• The leptons carry essentially all the available energy, and the nuclear recoil is 



• The transition involves the 0"*" ground state of the initial nucleus and, in almost all 
cases, the 0'*' ground state of the final nucleus. For some isotopes, it is energetically 
possible to have a transition to an excited O"*" or 2+ final state(^), even though they 
are suppressed because of the smaller phase space available; 

• Both processes are second-order weak processes, i.e. their rate is proportional to 
Gp, where Gf is the Fermi constant. They are therefore inherently slow. Phase 
space considerations alone would give preference to the /3/30nu mode which is, 
however, forbidden by total lepton number conservation. 

The distinct features are: 

• In the I3f52v mode the two neutrons undergoing the transition are uncorrelated (but 
decay simultaneously), while in the /3f30h' mode the two neutrons are correlated; 

• in the /3/32v mode, the sum electron kinetic energy Ti + T2 spectrum is continuous 
and peaked below (5/3/3/2, where Qpp is the Q-value of the reaction. In the fijiQv 
mode, since no light particles other than the electrons are emitted and given that 
nuclear recoil is negligible, the Ti + T2 spectrum is a mono-energetic line at Qpp, 
smeared only by the detector resolution. This is illustrated in fig. 7. 

In addition to the the two basic decay modes described above, several decay modes 
involving the emission of a light neutral boson, the Majoron (x"), have been proposed in 
extensions of the Standard Model, see sect. 3'4. 

While in the following we will focus on /3/?0i^ as defined in eq. (23), there are three 
closely related lepton number violating processes that can be investigated: 



Such processes are called double positron emission, single positron emission plus single 
electron capture (EC), and double electron capture, respectively. All three involve tran- 
sitions where the nuclear charge decreases (as opposed to increasing, as in (3/30i>) by two 
units. From the the theoretical point of view, the physics probed by /3+/3+Ot^, /3+ECOi^ 
and ECECOi^ is identical to the one probed by PfiQv. From the experimental point of 
view, however, P^Qv and /3+ECOi^ are less favorable than fijSQv because of the smaller 



(®) The transition to an excited O"*" final state has been observed for both ^"''Mo [49, 50, 51] 
and ^^°Nd [52]. 



negligible; 



(24) 
(25) 
(26) 



/3+ECOi/ 
ECECOi^ 



{Z,A) {Z -2,A) + 2 e+ 
e- + {Z,A)^ {Z -2,A) + e+ 
2 e- + {Z,A) {Z-2,A)* 
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(Ti+T2)/Q 



Fig. 7. - Spectra for the sum kinetic energy Ti + T2 of the two electrons, for different /3/3 modes: 
(3p2i/, PpOv, and /3/3 decay with Majoron emission. 



phase space available. On the other hand, the process ECECOi^ is gaining some atten- 
tion recently as a promising (but still much less developed) alternative to /3/30z^, since a 
resonant enhancement of its rate can in principle occur [53] . 

In the following, the neutrinoless mode (SpOv is discussed in more detail, from both 
the theoretical and experimental point of views. 

3'2. The black box theorem . ~ In general, in theories beyond the Standard Model 
there may be several sources of total lepton number violation which can lead to PjdQv. 
Nevertheless, as it was first pointed out in reference [54], irrespective of the mechanism, 
PPOv necessarily implies Majorana neutrinos. This is called the black box (or Schechter- 
Valle) theorem. The reason is that any AL 7^ diagram contributing to the decay would 
also contribute to the (e, e) entry of the Majorana neutrino mass matrix, (my)ee- This 
is shown in fig. 8, where a, De — i^e transition, that is a non-zero {m^)ee, is induced as a 
consequence of any AL 7^ operator responsible for /3f30i/. 

From a quantitative point of view, however, the diagram in fig. 8 corresponds to a tiny 
mass generated at the four-loop level, and is far too small to explain the neutrino mass 
splittings observed in neutrino oscillation experiments [55]. Other, unknown, Majorana 
and/or Dirac mass contributions must exist. As a consequence, therefore, the black box 
theorem says nothing about the physics mechanism dominating a jSpOv rate that is large 
enough to be observable. The dominant mechanism leading to j3l3Qv could then either be 
directly connected to neutrino oscillations phenomenology, or only indirectly connected 
or not connected at all to it [56]. The former case is realized in the standard P^Qv 
mechanism of light neutrino exchange, discussed in sect. 3'3. The latter case involves 
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Fig. 8. - Diagram showing how any neutrinoless double beta decay process induces a iz-to-i/ 
transition, that is, an effective Majorana mass term. This is the so-called black box theorem 
[54]. 



alternative /SpOiy mechanisms, briefly outlined in sect. 3'4. 

3'3. The standard Pj3Qv mechanism: light Majorana neutrino exchange. - Neutrinoless 
double beta decay can arise from a diagram (fig. 9) in which the parent nucleus emits a 
pair of virtual W bosons, and then these W exchange a Majorana neutrino to produce 
the outgoing electrons. The rate is non-zero only for massive, Majorana neutrinos. The 
reason is that the exchanged neutrino in fig. 9 can be seen as emitted (in association 
with an electron) with almost total positive helicity. Only its small, 0{m/E), negative 
helicity component is absorbed at the other vertex by the Standard Model clcctrowcak 
current. Considering that the amplitude is in this case a sum over the contributions of 
the three light neutrino mass states i^i , and that is also proportional to C/J^ , we conclude 
that the modulus of the amplitude for the Pj3Qi' process must be proportional in this 
case to the effective neutrino Majorana mass: 



(27) 



In other words, the effective neutrino Majorana mass corresponds to the modulus of the 
(e, e) element of the neutrino mass matrix of eq. (21), mjjp = | {m^)^^ \ . 

In the case where light Majorana neutrino exchange is the dominant contribution to 
/3p0iy, the inverse of the half-life for the process can be written as [57]: 

(28) -^=G'-^{Q,Z)\M''^\'ml^, 

^1/2 

where G^'^{Q, Z) is a phase space factor that depends on the transition Q- value and on 
the nuclear charge Z, and M^'^ is the nuclear matrix element (NME). The phase space 
factor can be calculated analytically, in principle with reasonable accuracy (^). The NME 
is evaluated using nuclear models, although with considerable uncertainty (see sect. 4). 



ij) An accurate description of the effect of the nuclear Coulomb field on the decay electron 
wave-functions is, however, required. 
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Fig. 9. - The standard mechanism for PPOu decay, based on light Majorana neutrino exchange. 



In other words, the value of the effective neutrino Majorana mass m^p in eq. (27) can 
be inferred from a non-zero Pj3Qv rate measurement, albeit with some nuclear physics 
uncertainties. Conversely, if a given experiment does not observe the PPQv process, the 
result can be interpreted in terms of an upper bound on mpp. 

If light Majorana neutrino exchange is the dominant mechanism for PPOv, it is clear 
from cq. (27) that f3/30i' is in this case directly connected to neutrino oscillations phe- 
nomenology, and that it also provides direct information on the absolute neutrino mass 
scale, as cosmology and (3 decay experiments do (see sect. 2'1). The relationship between 
nipp and the actual neutrino masses is affected by: 

1. the uncertainties in the measured oscillation parameters; 

2. the unknown neutrino mass ordering (normal or inverted); 

3. the unknown phases in the neutrino mixing matrix (both Dirac and Majorana). 

For example, the relationship between mpp and the lightest neutrino mass mught 
(which is equal to mi or in the normal and inverted mass ordering cases, respectively) 
is illustrated in fig. 10. This graphical representation was first proposed in [58]. The width 
of the two bands is due to items 1 and 3 above, where the uncertainties in the measured 
oscillation parameters (item 1) are taken as Scr ranges from a recent global oscillation fit 
[3]. Figure 10 also shows an upper bound on miight trom cosmology (mught < 0.43 eV), 
also shown in fig. 2, and an upper bound on m^g^ from current I3f30v data (m^^ < 
0.32 eV), which we will discuss in sect. 3'5. As can be seen from fig. 10, current /3^0i/ 
data provide a constraint on the absolute mass scale mught that is almost as competitive 
as the cosmological one. 

In figs. 2 and 10, we have shown only upper bounds on various neutrino mass combi- 
nations, coming from current data. The detection of positive results for absolute neutrino 
mass scale observablcs would open up the possibility to further explore neutrino prop- 
erties and lepton number violating processes. We give three examples in the following. 
First, the successful determination of both mp in cq. (3) and mpp in eq. (27) via /3 
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flight (eV) 

Fig. 10. - The effective neutrino Majorana mass mpp as a function of tlie liglitest neutrino 
mass, miight- The red (green) band corresponds to the normal (inverted) ordering, respectively, 
in which case miight is equal to mi (ma). The vertically-excluded region comes from cosmological 
bounds, the horizontally-excluded one from PPQv constraints. 



and PPOv decay experiments, respectively, can in principle be used to determine or con- 
strain the phases Ui [59]. Second, measurements of mp or mcosmo in eq. (2) may yield 
a constraint on miight that is inconsistent with a mp^ upper limit. In this case, the 
non-observation of PfiOv would suggest that neutrinos are Dirac particles. Third, mea- 
surements of mp or TOcosmo may yield a constraint on miight that is inconsistent with a 
measured non-zero mpjs. This scenario would demonstrate that additional lepton number 
violating physics, other than light Majorana neutrino exchange, is at play in the /3/30iy 
process. We briefly describe some of these possible PfSOv alternative mechanisms in the 
following. 

3'4. Alternative f3/30i> mechanisms. - A number of alternative P/30v mechanisms have 
been proposed. For an excellent and complete discussion of those, we refer the reader to 
[56]. The realization of (3f5Qu can differ from the standard mechanism in one or several 
aspects: 

• The Lorcntz structure of the currents. Positive chirality currents mediated by a Wr 
boson can arise, for example, in left-right symmetric theories. A possible diagram 
involving positive chirality current interactions of heavy Majorana neutrinos Ni is 
shown in fig. 11(a). 

• The mass scale of the exchanged virtual particles. One example would be the 
presence of "sterile" (that is, described by positive chirality fields) neutrinos, cither 
light or heavy, in the neutrino propagator of fig. 9, in addition to the three light, 
active, neutrinos we are familiar with. Another example would be the exchange of 
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Fig. 11. - Examples of non-standard mechanism for PjBQv: (a) heavy neutrino exchange with 
positive chirality currents [60]; (b) neutralino exchange in R-parity violating supersymmetry 
[61]; (c) Majoron emission [62]. 



heavy supersymmetric particles, as in fig. 11(b). 

• The number of particles in the final state. A popular example involves decay modes 
where additional Majorons, that is very light or massless particles which can couple 
to neutrinos, are produced in association with the two electrons (see fig. 11(c)). 

In non-standard jS/SOi' mechanisms, the scale of the lepton number violating physics 
is often larger than the momentum transfer, in which case one speaks of short-range 
processes. This is in contrast to the standard (ipOv mechanism of light Majorana neutrino 
exchange, in which case the neutrino is very light compared to the energy scale, resulting 
in a long-range process. Non-standard and long-range jSpQu processes are, however, also 
possible. 

In general, several contributions to the total PfiQv amplitude can add coherently, al- 
lowing for interference effects. Neutrinolcss double beta decay obscrvables alone may be 
able to identify the dominant mechanism responsible for jSpOv. We give two examples. 
First, if Majorons are also emitted in association with the two electrons, energy conser- 
vation alone requires the electron kinetic energy sum Ti -f T2 to be a continuous spectrum 
with Qpp as endpoint. This spectrum is potentially distinguishable from the jipiv one 
(see fig. 7), provided that the Majoron-neutrino coupling constant is large enough. Sec- 
ond, if positive chirality current contributions dominate the PfiQv rate, electrons will be 
emitted predominantly as positive helicity states. As a consequence, both the energy 
and angular correlation of the two emitted electrons will be different from the ones of 
the standard P^Qv mechanism. A detector capable of reconstructing individual electron 
tracks may therefore be able to distinguish this type of non-standard f5f5Qv mechanism 
from light Majorana neutrino exchange (see, for example, [63]). 

3'5. Existing experimental results. - Neutrinolcss double beta decay searches have been 
carried out over more than half a century, exploiting the same experimental techniques 
used for measuring the two-neutrino mode rate (see also sect. 6'1). Several /3/3 emitting 
isotopes have been investigated, as shown in table III. 

The most sensitive limit to date was set by the Heidelberg-Moscow (HM) experiment 
[1]: T^J^(^^Gc) > 1.9 X 10^5 years (90% CL), corresponding to an effective Majorana 
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Isotope 


rO/2 (years) 


Experiment 


«Ca 


> 5.8 X 10^2 


ELEGANT [64] 




> 1.9 X 10^^ 


Heidelberg-Moscow [1] 




> 3.6 X 10^3 


NEMOS [65] 


86 Zr 


> 9.2 X 10^1 


NEM03 [38] 


Mo 


> 1.1 X lO^"* 


NEM03 [65] 




> 1.7 X 10^^ 


Solotvina [42] 


130 rpg 


> 2.8 X 10^* 


CUORICINO [66] 




> 4.5 X 10^^ 


DAMA [67] 




> 1.8 X 10^2 


NEM03 [46] 



Table III. - Current best limits on the half-life of PPOu processes for the most interesting 
isotopes. All values are at 90% C'L. 



mass bound of m^/j < 0.32 cV. A subset of this collaboration claimed to observe 
evidence for a f3(30i/ signal, with a best value for the half-life of 1.5 x 10^^ years [68]. 
The claim is very controversial [69]. Also, a subsequent re-analysis by the same group 
updated this result to (2.23^0 3^) x 10^^ years [70], resulting in an effective Majorana 
mass of about 0.30 eV. The authors claim a statistical significance for the evidence of 
6(7, and do not present any systematic uncertainty analysis. It should be mentioned 
that a different germanium detector, IGEX. did not observe any evidence for /3/30i^, and 
reported a lower boimd on the half-life of: T°Y^C^Ge) > 1.57 x 10^^ years (90% CL) [71]. 



4. Calculating nuclear matrix elements 

All nuclear structure effects in /3(30i/ are included in the nuclear matrix element 
(NME). Its knowledge is essential in order to relate the measured half-life to the neutrino 
masses, and therefore to compare the sensitivity and results of different experiments, as 
well as to predict which are the most favorable nuclides for /SpOv searches. Unfortunately, 
NMEs cannot be separately measured, and must be evaluated theoretically. 

In the last few years the reliability of the calculations has greatly improved, with sev- 
eral techniques being used, namely: the Interacting Shell Model (ISM) [72]; the Quasi- 
particle Random Phase Approximation (QRPA) [48]; the Interacting Boson Model (IBM) 
[73]; and the Generating Coordinate Method (GCM) [74]. Before briefly reviewing the 
different approaches, we discuss the ingredients that are common to all. 

It is beyond the scope of this review, and beyond our expertise, to provide a com- 
plete derivation of NME calculations. Here, we limit ourselves to outline the theoretical 
framework used to carry out the calculations, the approximations used, and the most sig- 
nificant differences among the different approaches. For more details, the reader should 
refer to [75, 59, 48], where most of the material covered below was taken from. 

4'1. Common elements in calculations. - The rate for the PPOv process can be written 

as: 



(29) 



/2\ 



E 

spins " 



\Zo,\^SiE,i+E,2 + Ef-AU) 



27r3 27r3 
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In this formula, i?i(2) and Pi{2) £^re the total energies and momenta of the electrons, 
Ef{Mi) is the energy of the final (mass of the initial) nuclear state, and Zq^ is the 
reaction amplitude, to be evaluated in time-dependent perturbation theory to second 
order in the weak interaction (that is, to second order in the Fermi constant Gp)- The 
reaction amplitude can be factorized into the product of a leptonic and a hadronic part. 
Assuming that the decay is mediated solely by the exchange of light neutrinos (standard 
PPQv mechanism, see sect. 3'3), the leptonic part can be written as the product of two 
negative chirality currents. After substitution for the neutrino propagator, the lepton 
amplitude acquires the form: 

(30) / E (0e--(-^)e(x)7,(l - 75)^^?^(1 - l,)l.e^{y) Uj, , 

where the integral is over the 4-momentum transfer q (that is, the momentum of the 
virtual neutrino), the sum k is over the three neutrino mass eigenstates of mass m/c, the 
mixing matrix elements Uek specify the electron flavor content of the mass states, and 
e{x) and e'^{y) are the electron creation operators. This lepton part implies a contraction 
over the two neutrino operators, which is allowed only if neutrinos are Majorana particles. 
Furthermore, from the commutation properties of the gamma matrices it follows that 

(31) ^(1 - l^Kq'lp + mk){l - 75) = ^(1 - 75) 

From eqs. (30) and (31), we obtain that the decay amplitude for purely negative chirality 
lepton currents is proportional to the neutrino Majorana mass U'^^jrik, as discussed 
in the previous section. 

Integration over the virtual neutrino energy leads to the replacement of the propagator 
(g^ — m\)~^ by the residue tt/w^ with = \/ + m| . Integration over the space part 
dq leads to an expression representing the effect of the neutrino propagation between the 
two nucleons. This expression has the form of a neutrino potential, H[r), where r < R, 
and R is the nuclear radius, R = 1.2 A^^^ fm. H{r) appears in the corresponding nuclear 
matrix elements, introducing a dependence of the transition operator on the coordinates 
of the two nucleons, as well as a weak dependence on the excitation energy Em — Ei of 
the virtual state in the odd-odd intermediate nucleus. 

The momentum of the virtual neutrino is determined by the uncertainty relation 
g ~ l/r. Here r is a typical spacing between two nucleons, r ~ 2 — 3 fm. Therefore the 
momentum transfer is q ^ 100 — 200 MeV. For light neutrinos the neutrino mass ruj 
can then be safely neglected in the potential H{r). Also, given the large value of q, the 
dependence on the difference of nuclear energies Em — Ei is weak. One can then neglect 
the variation of the energy from state to state when integrating over the virtual neutrino 
energies. In this closure approximation the contributions of the two electrons are added 
coherently, and the neutrino potential is of the form: 

(32) H{r) = -$(wr) 

r 

The nuclear radius R in eq. (32) is introduced in order to make the potential H dimen- 
sionless. A corresponding 1/R^ factor compensates for this auxiliary quantity in the 
phase space formula (see eq. (40) below). Also, in eq. (32), $(ajr) < 1 is a relatively 
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slowly varying function of r. A typical value of H{r) is larger than unity, but less than 
5-10. 

To go any further, we need an expression for the hadronic current. In the impulse 
approximation, the hadronic current is obtained from that of free nucleons. The latter 
can be written as 

(33) J^t = *T+ [gviq^h" - 9A{q^)Yl5 - gp{q^)q''l5] * , 

where nip is the nucleon mass, ^ is a nucleon field, r+ = (ti + iT2)/2, ri(2) are the 
isospin Pauli matrices, and gv, gA and gp are the so-called vector, axial-vector and 
induced pseudoscalar form factors, parametrizing the composite structure of nucleons. 
Since in PPOv decay we have q^ ^ qQ, we take ~ —q^- 

The q^ dependence of the vector and axial- vector form factors is parametrized via the 
usual dipole approximation 

(34) 5v(g')-gy/(l + 'f'A4)', 9Aiq^)^9A/il + qVMlf , 

with gv ~l,gA^ 1-25, My ~ 850 MeV, and Ma ^ 1090 MeV. It is also customary to 
use the Goldberger-Treiman relation for the induced pseudoscalar term 

(35) gp{q^) = 2mpgA{q^)/{q'- + ml) . 

where m.^ is the pion mass. For the ground state to ground state transitions, i.e. — > 
0^, it is sufficient to consider s-wavc outgoing electrons {long-wave approximation), and 
the nonrelativisitic approximation for the nucleons. The rate then takes the form given 
in eq. (28), and repeated here for convenience: 

(36) [^J^'^G^^iQ^Z) lA/O-'l' ml, , 

where Q = Mi — Ef, G^^{Q,Z) comes from the phase-space integral, and the nuclear 
matrix element turns into a sum of the Gamow- Teller and Fermi nuclear matrix elements, 
where : 

(37) M"''^(^)'(^M°^^- Ja/O^) 
with, to first order: 

(38) M^^ = {f\Y,H{r)TtT^\i) 

a,b 

and 

(39) M^Jt = (/I E ^(^)^'' • '^'^^aT^ K) • 

a,b 

In these equations the neutrino potential H{r) is of the form defined in eq. (32) (for 
explicit realizations of the neutrino potential see, for example, [48]), (7a{b) are spin Pauli 
matrices, and |/) {\i)) are the final (initial) nuclear states. In contrast to Ii(i2v, which 
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Fig. 12. - The phase space factor for all /3/3 emitters with Q > 2 MeV. Values taken from [76, 77] 



involves only Gamov- Teller transitions through intermediate 1+ states (because of low 
momentum transfer), the nuclear matrix element for (3j3Qv involves all multipolarities 
in the intermediate odd-odd (A, Z + I) nucleus, and contains both a Fermi (F) and a 
Gamov-Teller (GT) part. 

In eq. (36) the explicit form of the phase-space integral is: 

G-(Q,Z).(G.W(1)^^^. 

(40) (^J F{Z,Eel)F{Z,Ee2)PelPe2EelEe2S{Eo-Eel-Ee2)dEeldEe2^ ■ 

where Eq = Q + 2mp = Mi — Mf is the available energy, pei (Pe2) a-re the electron 
3-momenta, F{Z, E) is the Fermi function that describes the nuclear Coulomb effect on 
the outgoing electrons, and Z is the charge of the daughter nucleus. 

If an accurate result is required, the relativistic form of the function F{Z, E) must 
be used and a numerical evaluation is necessary [57]. The phase space factor for all 
PP emitters with Q > 2 MeV are given in fig. 12. 

For a qualitative picture, one can use the simplified nonrelativistic Coulomb expres- 
sion, the so-called Primakoff- Rosen approximation [78]: 



(41) F{Z,E)^^-^^ 
In this approximation, G^'' is independent of Z: 

(42) + ^ + 
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r [fm] 

Fig. 13. - The dependence of M^" as a function of tlie distance r among the two neutrons 
participating in P/SOu, in ^®Ge. The four curves show the effects of different treatments of 
nucleon-nucleon short-range correlations [48]. 



where Eq is expressed in units of electron mass. Notice that the phase space dependence 
of the /3/3Qi' mode goes with Eq, while the phase space of the corresponding two- neutrino 
mode goes with Eq^. That is, based on phase space considerations alone, the /3/30j^ mode 
would be much faster than the /3/32:/ mode, if the neutrino mass were of the order of the 
electron mass. 

In addition to the total NME in eq. (37), the different nuclear structure approaches 
discussed in sect. 4'2 allow to predict also what are the typical distances among the two 
decaying nucleons that contribute the most to Af°''. The result is shown in fig. 13. One 
can see that only relatively short distances, r < 2 — 3 fm, contribute significantly. In 
other words, essentially only the nearest neighbour neutrons undergo P^Qv transitions. 
Although these short distances justify the above-mentioned closure approximation, the 
fact that two nucleons strongly repel each other for distances r < 0.5 — 1.0 fm should 
be taken into account. The nuclear wavefunctions computed according to the methods 
described in sect. 4'2 do not take such effect into account. The usual and simplest 
way to include this effect is by introducing a plienomenological function in eq. (37). A 
popular procedure to obtain such function is based on the Unitary Correlation Operator 
Method (UCOM) [79]. This procedure reduces the value of M^" in fig. 13 by only about 
5%. However, other prescriptions for short-range correlations introduce a much more 
significant (of order 20-25%) reduction in M'^" . 

4'2. The different nuclear structure approaches. - In order to compute the l3/30i> 
decay rates for a given neutrino mass, we need to evaluate the initial and final state 
wavefunctions \i) and )/), and the nuclear matrix elements connecting the two in eq. (37). 
Given the complicated nuclear many-body nature of the problem, this calculation cannot 
be done exactly, and some approximations need to be introduced. Different nuclear 
physics approaches have been used to this end. Only a very schematic description of 
those will follow. 
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Fig. 14. - Low-lying energy levels in a single-particle shell model with an oscillator potential 
(with a small negative P term) without spin-orbit (left) and with spin-orbit (right) interaction. 
The number to the right of a level indicates its degeneracy, (2j-f 1). The boxed integers indicate 
the magic numbers. 



4'2.1. The Interacting Shell Model. In the Interacting Shell Model (ISM) [80], aU 
microscopic calculations are based on the Independent Particle (Shell) Model (IPM). 
The basic premise of such a model is that the nucleons are moving independently in a 
mean field with a strongly attractive spin-orbit term: 

(43) U{r) = ^hujr^ +Dl^ + Cl- s. 

where the harmonic oscillator (plus the surface correction D I part describes the 
bound nucleon nature of the problem, and the spin-orbit part C I ■ s is added to give the 
proper separation of the subshells and explain the nuclear magic numbers, i.e. specific 
values of the number of protons Z and neutrons N {N or Z = 2, 8, 20, 28, 50, 82, 126) 
accounting for the existence of shell closures at those occupation numbers. In Eq. 43, I 
is the orbital angular momentum, and s the spin of single nucleons. The effect of the 
spin-orbit potential on the nuclear energy levels is schematically shown in Fig. 14. 

As the number of protons and neutrons depart from the magic numbers, it becomes 
indispensable to include the "residual" two-body nucleon interaction among nucleons. 
This point marks the passage from the IPM to the ISM model. This residual interaction 
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contains both a kinetic (/C) and a potential (V) term: 



(44) 




i<j k<i 



that adds one or two particles in orbits of total angular momentum i, j and removes one 
or two from orbits k, I, subject to the Pauli principle ({a|aj} = Sij). While complicated 
in practice, this approach is conceptually simple: given a good enough residual interac- 
tion Vijki, the problem is reduced to diagonalizing a matrix in a sufficiently large basis 
("valence space"). In this framework, a limited valence space is used but all configura- 
tions of valence nucleons are included. The ISM describes well properties of low-lying 
nuclear states. 

4'2.2. The Quasiparticle Random Phase Approximation. The basic idea behind the 
proton-neutron Quasiparticle Random Phase Approximation (QRPA) [48] is that the 
most important part of the residual interaction among nucleons is the pairing force. The 
pairing force accounts for the tendency of nucleons to couple pairwise to especially stable 
configurations, i.e. into nuclei with even TV, even Z. This force favors the coupling of 
neutrons with neutrons, and protons with protons, so that the orbital angular momentum 
and spin of each couple adds to zero. As the result of the pairing force, the nuclear ground 
state is mainly composed of Coopcr-like pairs of neutrons and protons coupled to = 0'*" 
total angular momentum. In QRPA, the nucleon pairing is introduced via the BCS theory 
of superconductivity. A unitary (Bogoliubov) transformation is first performed to change 
from a particle to a quasiparticle basis. Quasiparticles are generalized fermions which are 
partly particles (with probability , where j is the single-particle orbital the quasiparticle 
belongs to, see fig. 15) and partly holes (with probability v'j). Quasiparticles are just a 
mathematical construct to account for pairing between like nucleons in a simple fashion 
while retaining the simplicity of the independent particle model, since the quasiparticles 
are kept, to first order, independent. This transformation smears out the nuclear Fermi 
surface over several orbitals, for both protons and neutrons, as shown schematically in 



Once the problem has been transformed into the simpler quasiparticle basis, the 
QRPA goal is to evaluate the transition amplitudes associated with charge changing one- 
body operator T'^^^ connecting the 0^ vacuum of quasiparticles in the even-even nucleus 
with any of the J'^ excited states in the neighboring odd-odd nuclei. Such states are 
described as harmonic oscillations above this vacuum. The creation of such particle-hole 
pairs (or phonons) from a BCS-only vacuum would, however, overestimate this transition 
amplitude. The QRPA is the simplest theory which admits the possibility that the ground 
state is not of purely independent quasiparticle character, but may contain correlations. 
As a consequence, two-particle, two-hole excitations are included in the QRPA vacuum 
state, as opposed to the BCS vacuum. The transition amplitude is then modified as 
needed, since the creation of a particle-hole pair from the BCS vacuum (the so-called 
forward- going amplitude X) can lead to the same final state J'^ as the destruction of a 
particle-hole pair from a two-particle, two-hole excitation (the backward- going amplitude 
Y). The amplitudes X and Y as well as the corresponding energy eigenvalues Wm are 
determined by solving the QRPA equations of motion for each J'^: 



fig. 15. 



(45) 
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(a) 



(b) 



Fig. 15. - Occupation probabilities of single particle orbitals in (a) the Independent Particle 
Model, and (b) with the inclusion of pairing forces. The arrows indicate possible excitations of 
the nucleus induced by transfer of a particle from a (partially) occupied orbital to a (partially) 
unoccupied orbital. Adapted from [81]. 



In eq. (45) the terms A and B depend on the interaction matrix elements between 
quasi-particle configurations. They can be written in terms of particlc-holc (p-h) and 
particle-particle (p-p) matrix elements. Customarily, these interaction matrix elements 
are multiplied by adjustable coupling constants gph and gpp, respectively. If a realistic 
nucleon-nucleon interaction is used, then the values of these constants are gph — Qpp — 1. 
These particle-particle interactions enhance the backward-going amplitude Y , thereby 
reducing the transition amplitude. 

4'2.3. The Generating Coordinate Method. A nucleon coupling scheme that competes 
with nucleon pairing (see sect. 4'2.2) to fix the equilibrium shape of a nucleus, and its 
collective motion, is the so-called aligned coupling scheme. In this scheme, each nucleon 
has the tendency to align its orbit with the average field produced by all other nucleons. 
This preferentially gives rise to nuclei with deformed equilibrium shapes and collective 
rotational motion. A common representation of the shape of these nuclei is that of an 
ellipsoid. The quadrupole deformation parameter j3 is related to the eccentricity of the 
ellipse: /? 7^ represents a non-spherical nucleus, with /? > (/3 < 0) corresponding to a 
prolate (oblate) ellipsoid. 

Nuclear collective rotors arc associated with "intrinsic states" very well approximated 
by deformed mean field determinants, that is antisymmetrized products of independent 
particle wavefunctions. Nuclear wavefunctions of this type represent the basic assumption 
of the Hartree-Fock self-consistent field theory (see, for example, [81]), an approximation 
for reducing the problem of many interacting particles to one of non-interacting particles 
in a mean field. 

In the Generating Coordinate Method with Particle Number and Angular Momentum 
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Projected product-type wave functions (GCM-PNAMP, or GCM in short) [82], one starts 
by building a set of Hartree-Fock-Bogoliubov (HFB) intrinsic axial symmetric wave func- 
tions |(/)/3) along the quadrupole deformation parameter /3. These HFB intrinsic states are 
found by solving the so-called constrained Particle Number Variation After Projection 
equations (PN-VAP), ^(£'^■•^(10^))) = 0. This is a variational equation constrained to a 
fixed value of the quadrupole deformation /3. The Gogny DIS interaction is used as the 
underlying nucleon-nucleon interaction. Exact eigenstates can be found by projecting 
from the HFB wavefunctions the components of well-defined angular momentum, proton 
number and neutron number. The initial and final ground states |0^) can therefore be 
written as GCM wavefunctions: 



where P^^^, P^ and P^ are the angular momentum (/ for axial symmetric wave- 
functions), neutron number and proton number projection operators. In this method, 
one starts with the projected HFB approach, but allows for admixtures of different de- 
formations P, as described by eq. (46). The coefficients of this admixture arc found 
by solving the so-called Hill- Wheeler- Griffin (HWG) equation of generator coordinates 
(see, for example, [81]). 

4'2.4. The Interacting Boson Model. A somewhat intermediate path between the 
"microscopic" view of nuclear structure (ISM) and the "collective" views (QRPA, GCM) 
mentioned above was opened by the Interacting Boson Model, IBM [83]. In the inter- 
acting boson model, collective excitations of nuclei are described by bosons, and the 
microscopic foundation of such collective nuclear states is rooted in the shell model. 
As the number of valence nucleons increases, the direct application of the shell model 
becomes prohibitively difficult, and some approximation is needed. First, one usually 
assumes that the closed shells are inert. Second, one assumes that the important par- 
ticle configurations in even-even nuclei are those in which identical particles are paired 
together in states with total angular momentum and parity = 0"*" or = 2"*". Third, 
one treats the pairs as bosons, much in the same way as Cooper pairs in a gas of electrons. 

If one retains all three approximations, one is led to consider a system of interacting 
bosons of two types, proton bosons and neutron bosons. The proton (neutron) bosons 
with = 0+ are denoted by (sjy), the ones with = 2+ are denoted by {di,). 
The multitude of shells which appears in the shell model is then reduced to the simple 
s-shell (J = 0) and the d-shell (J = 2). The number of proton (A'^) and neutron [Nif) 
bosons is counted from the nearest closed shell, i. e. if less (more) than half of the shell 
is full, A^^(i/) is taken as the number of particle (hole) pairs. 

All fermionic operators, for example the operators yielding the Gamow- Teller and 
Fermi nuclear matrix elements in eqs. 38 and 39, are similarly mapped into bosonic 
operators by the Otsuka, Arima and lachello (OAI) method [84]. Using this method one 
is assured that the matrix elements between fermionic states in the collective subspace 
are identical to the matrix elements in the bosonic space [73]. 

A realistic set of wavefunctions of even-even nuclei with mass A > 60 is provided 
by the proton-neutron IBM-2 [83]. The wavefunctions arc generated by diagonalizing 
the IBM-2 Hamiltonian. These IBM-2 wavefunctions provide an accurate description 
of many properties (energies, electromagnetic transition rates, quadrupole and magnetic 
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Fig. 16. - Recent NME calculations from different techniques (GCM [74], IBM [73], ISM [86, 87], 
QRPA(J) [88], QRPA(T) [89, 90, 91]) with UCOM short range correlations. All the calculations 
use gA = 1.25; the IBM-2 results are multiplied by 1.18 to account for the difference between 
Jastrow and UCOM, and the RQRPA are multiplied by 1.1/1.2 so as to line them up with 
the others in their choice of ro = 1.2 fm. The shaded intervals correspond to the proposed 
physics- motivated ranges (see text for discussion). 



moments, etc.) of the final and initial nuclei. Using these wavefunctions, and the bosonic 
operators of the OAI method, it is possible to calculate /3/30i^ NMEs (for details, see [73]). 

4'3. Quantifying uncertainties in NME calculations. - Figure 16 shows the results of 
the most recent NME calculations with the methods described in sect. 4'2. We can see 
that in most cases the results of the ISM calculations are the smallest ones, while the 
largest ones may come from the IBM, QRPA or GCM. For a detailed study quantifying 
the spread of NME results resulting from different calculations, we refer the reader to 
ref. [85]. Shall the differences between the different methods in sect. 4'2 be treated as 
an uncertainty in sensitivity calculations? Should we assign an error bar to the distance 
between the maximum and the minimum values? This approach, we argue, does not 
reflect the recent progress in the theoretical understanding of the treatment of nuclear 
matrix elements. In quantifying the uncertainties in NME calculations, we follow [92]. 

Each one of the major methods has some advantages and drawbacks, whose effect 
in the values of the NME can be sometimes explored. The clear advantage of the ISM 
calculations is their full treatment of the nuclear correlations, while their drawback is that 
they may underestimate the NMEs due to the limited number of orbits in the affordable 
valence spaces. It has been estimated [93] that the effect can be of the order of 25%. 
On the contrary, the QRPA variants, the GCM in its present form, and the IBM are 
bound to underestimate the multipole correlations in one or another way. As it is well 
established that these correlations tend to diminish the NMEs, these methods should 
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tend to overestimate them [72, 94]. 

With these considerations in mind, physics-motivated ranges (PMR) of theoretical 
values for ^^Ge, ^^^Se, i^^Te, "^Xe and ^^^Nd NMEs have been proposed in [92]. In 
quantifying the uncertainties, the results of the major nuclear structure approaches which 
share the following common ingredients were considered: (a) nucleon form factors of 
dipole shape, see eq. (34); (b) soft short-range correlations computed with the UCOM 
method; (c) unquenched axial coupling constant gA = 1.25; (d) higher order corrections 
to the nuclear current [77] accounted for; and (e) nuclear radius R = tq A^^^, with 
To = 1.2 fm [95]. Therefore, the remaining discrepancies between the diverse approaches 
arc solely due to the different nuclear wavefunctions that they employ. The uncertainties 
in NME calculations for '^^Ge, ^^Se, ^^^Te, ""^Xe and ^^^Nd are shown as grey bands in 
fig. 16, and are in the 20-30% range. 



5. Ingredients for the ultimate P/SOiy experiment 

The discovery of l3j3Qv would represent a substantial breakthrough in particle physics. 
A single, unequivocal observation of the decay would prove the Majorana nature of 
neutrinos and the violation of lepton number. Alas, that is not, by any means, an easy 
task. The design of a detector capable of identifying efficiently and unambiguously such 
a rare signal represents a major experimental problem. 

To start with, one needs a large mass of the scarce /3/3 isotope in order to probe in 
a reasonable time the extremely long lifetimes expected. For instance, for a Majorana 
neutrino mass of 50 meV, it can be estimated using eq. (28) and a sound assumption 
for the NMEs that half-lives in the range of 10^^-10^'' years must be explored (i.e., 17 
orders of magnitude longer than the age of the universe!). A better sense of what such 
extremely long half-lives mean can be grasped with a simple calculation. Consider the 
radioactive decay law in the approximation T1/2 ^ i, where t is the exposure time; in 
that case, the expected number of events is given by 



(47) A^,,o. = log2 ^''^^•^^ 



Wpp TO;,' 

where M^^ is the mass of the emitting isotope, Na is the Avogadro constant, Wpjj 
is the molar mass of the /3/3 isotope, and e is the signal detection efficiency. 

It follows from eq. (47) that, in order to observe (assuming perfect detection efficiency 
and no disturbing background) as little as one decay per year and assuming a Majorana 
neutrino mass of 50 meV (rjy'j ^ 10^^-10^'' years), "macroscopic" masses of /3/3 isotope 
of the order of 100 kg are needed. 

The situation becomes even more desperate when considering real experimental condi- 
tions. The background processes that can mimic a fipOv signal in a detector are copious. 
In the first place, the experiments have to deal with an intrinsic background, the l3/32v, 
that can only be distinguished by measuring the energy of the emitted electrons, since the 
neutrinos escape the detector undetected (see fig. 7). Good energy resolution is therefore 
essential to prevent the P/32v spectrum tail from spreading over the /3/30i^ peak. Never- 
theless, this energy signature could not be enough per se: a continuous spectrum arising 
from natural radioactivity can easily overwhelm the signal peak. Other signatures, like 
particle identification or the observation of the daughter nucleus, are a bonus to provide 
a robust result. 
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Several other factors such as detection efficiency or the scalabiHty to large masses 
must be taken into account as well when choosing the experimental technique. The 
simultaneous optimization of all these parameters is most of the time conflicting, and 
consequently, many different experimental approaches have been proposed and are under 
development. In order to compare their merits, a figure of merit, the experimental 
sensitivity to is normally used. We describe it below, followed by a discussion on 

the main parameters entering this figure. 

Before doing so, we note that our notation in eq. (47) — and in the rest of this review 
— differs from the usually adopted one, derived from source-equals-detector experimental 
configurations. In the source-equals-detector notation, one refers to the total active mass 
M of the detector, which is related to the mass Mpp in the /?/3 isotope via the following 
relationship: 



where W is the molecular weight of the molecule of the active material, a is the isotopic 
abundance of the candidate PjSQv nuclide, and rj is the number of (SpOv element nuclei 
per molecule of the active mass. For example, Te02 bolometric detectors with a natural 
isotopic abundance in ^-^"Tc are characterized by Wpp = 129.9 g/mol, W = 159.6 g/mol, 
a = 0.34167 and rj = 1, such that Mpf, = 0.278Af .(*) {^) 

5'1. Sensitivity of a PjSQv experiment. - All jSpQi/ experiments have to deal with 
non-negligible backgrounds, an only partially efficient PjSQv event selection, and more or 
less difficulties to extrapolate their detection technique to large masses. It is instructive, 
however, to imagine an ideal, background-free, experiment. If such an experiment, after 
running for an exposure Mpp ■ t, observes no events, it would report an upper limit in 
the l3j3Qv decay rate {T^J^)~^ , or possibly in the more relevant physical parameter mpp: 



where Ki is a constant that depends only on the isotope type, and on the details of the 
statistical method (and the confidence level) chosen to report such limit. Equation (49) 
follows directly from eqs. (28) and (47), see [92] for details. 

Let us now consider the sensitivity in the case of an experiment with background. In 
the large background approximation, the sensitivity as a function of the background rate 
h follows the classical limit: S{b) oc where b is the mean predicted background level. 



(*) To stress this somewhat unconventional mass notation and to avoid any confusion, we will 
make use in the following of kg/j^ as the mass unit to indicate one kilogram of /3/3 emitter mass. 
(®) As pointed out to us, in principle the best quantity to express the PPQv exposure, and 
the background rate per unit exposure and unit energy discussed below, is neither kg-year nor 
kg^^-year, but rather upp-yeax, where npp = Mpp ■ Na/Wpp is the number of moles of the 
nuclide. The reason is that 1 kg^/3 of, say, ^^Ge contains almost twice as many /3/3 nuclides as 1 
kgflfl of "°Nd. To avoid an even more "radical" departure from commonly employed units, we 
stick to kg^;3-year units in the following. 



(48) 



(49) 
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In this limit, the to^^ sensitivity can be written as 



(50) 



61/2 



e ■ Mpp ■ t 



where K2 is a constant depending on the isotope. If the background h is proportional to 
the exposure M^p ■ t and to an energy window /S.E around Qp^'. 



In short, the background limits dramatically the sensitivity of a double beta decay 
experiment, improving only as [Mpp ■ t)^^/^ instead of the [Mjip ■ t)^^^^ expected in the 
background-free case. 

Two aspects of cq. (51), and in particular of our definition of the background rate c, 
deserve further clarification. First, for a given background level 6, the background rate 
c will in general depend on the choice of the energy window AE. This is the case if the 
background energy spectrum around QjSfs is not flat. In the following, all background rate 
values refer to a AE choice of 1 FWHM energy resolution total width, that is computed 
for background events whose reconstructed energy falls in the [Qpfi — 0.5 • FWHM, + 
0.5 • FWHM] range. Similarly, the background rate c will in general depend on the mass 
Mjjp of the /3/3 emitting material considered. This is the case for backgrounds that are not 
uniformly distributed within the active mass, such as surface contaminations of materials 
or backgrounds that are of external origin. As already assumed in deriving eq. (52), all 
background rate values are relative to the total mass Mjsp appearing in the signal count 
rate computation of eq. (47). 

5'2. Choice of the (3/3 isotope. - In nature, 35 naturally-occurring isotopes are /3/3 
emitters. Which ones are the most favorable in terms of f3f30v searches? 

Let us start with considerations about the most favorable /3/30i^ phase space factors 
and nuclear matrix elements. We are interested in the isotopes that provide the highest 
(3j3Qv rate for the same m^p mass, or, in other words, those that minimize the constants 
Ki and/or K2 appearing in eqs. (50) and (52), respectively. To a first approximation, the 
phase space factor G^''{Q,Z) appearing in eq. (40) varies as see eq. (42). Isotopes 
with large Q-values are therefore favored. For this reason, only isotopes with > 2 
MeV are usually considered for /3f30v searches. The 11 isotopes satisfying this criterion 
are shown in fig. 12. The isotopes with the most favorable phase space factors are ^^^Nd, 
^®Ca and ^^Zr. As far as the nuclear matrix elements are concerned, variations from 
one isotope to another arc significantly smaller than G^'^{Q, Z) variations, as can be seen 
from fig. 16. Considering the relevant product of the phase space factor times the nuclear 
matrix element squared for the most promising f3P isotopes, we find variations of about 
a factor of 2 in mp^ sensitivity depending on the isotope and for an ideal experiment, as 
can be seen in fig. 17. From this figure, and from phase space factor and nuclear matrix 
element considerations alone, wc would conclude that ®^Se. I'^^Tc and i^°Nd would be 



(51) 



b = c- Mpii - t- AE 



with the background rate c expressed in counts/ (keV • kg • year), then: 



(52) 
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Fig. 17. - Sensitivity of ideal experiments at 90% CL for different /3/3 isotopes. Since the yields 
are very similar, the sensitivities of *^Se, ^''"Te and ^^'^Nd overlap . From reference [92]. 



preferable than '^Ge. However, other factors enter in the isotope choice, as discussed in 
this section. 

Another advantage in choosing a /3/3 isotope with a high Qpp value relies in background 
control. As we will see later, backgrounds to /3(30i/ searches from natural radioactivity 
populate the energy region below ~3 MeV. The possibility to use an isotope with Qis^ 
above these background energies is therefore desirable. 

One is also typically interested in choosing an isotope with a relatively slow l3(32i' 
mode. As the energy resolution degrades, the experiments are affected by /3/32i/ back- 
grounds in a more or less pronounced way, depending on the isotope. This is true unless 
the energy resolution of the experiment is truly excellent, in which case even relatively fast 
(3/321/ modes do not constitute a serious background to 13(301' searches. This is illustrated 
in fig. 18. In this figure, the m/j^ sensitivity (computed according to the prescription de- 
scribed later, in sect. 6T1) at 90% CL is shown for ideal experiments using five different 
isotopes as a function of FWHM energy resolution. The experiments, each assumed to 
use 100 kg^^ of /?/? emitter mass and to run for five years, are ideal in the sense of hav- 
ing perfect /3/30;/ efficiency and of being affected only by /3(32i' backgrounds. As fig. 18 
illustrates, and as far as the (3f32i' background is concerned and for the same moderate 
energy resolution (say, 5-10% FWHM), -'^'^^Xc is to be preferred over ^^Sc and ^'^"Nd, 
thanks to its much longer (3(321/ half-life (see tab. H). For experiments featuring excel- 
lent energy resolution, say < 2% FWHM, all experiments would operate in a essentially 
background-free regime, for the assumed 500 kg- yr exposure. As wc have seen above, in 
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Fig. 18. - Sensitivity to at 90% CL as a function of FWHM energy resolution, for ideal 
experiments using five different isotopes, each with 100 kg^^ of emitter mass and 5 years 
of data-taking. The experiments are assumed to have perfect efficiency and to be affected 
only by I3j32v backgrounds. In practice, experiments using ^^Ge and ^'^"Te always feature an 
excellent energy resolution and are therefore not affected by I3j32v backgrounds, hence only the 
background-free sensitivity limit is shown in those cases, with an arrow. 



this regime an isotope such as -^^°Nd is to be preferred over, say, ''^Ge(^°). In practice, 
however, other backgrounds are always present, typically creating a continuum through 
the region of interest, and a better resolution improves the experimental sensitivity even 
in the < 2% FWHM energy resolution range. 

Another factor entering in the /3/3 isotope choice has to do with how well understood 
the nuclear physics for that isotope is. As we have seen in sect. 4, the calculation of 
nuclear matrix elements is a very complicated task. In sect. 4'3, we have made an attempt 
at quantifying the uncertainties in the NMEs for various isotopes. Our conclusion is 
that no magic isotope exists, and uncertainties in the 20-30% range (according to our 
evaluation) exist for the five isotopes we have considered, '^^Ge, ^^Se, ^^Oxe, ise^e and 

150Nd. 

5'3. Isotope mass. - As explained above, large masses of isotope are needed to 
explore the expected half-lives. The previous generation of double beta decay experiments 
used masses of the order of 10 kg. New-generation experiments will range from tens of 



(^'') The cases of ^^Ge and ^^°Te are shown in fig. 18 only in this background-free regime, given 
that in practice experiments using such isotopes always feature excellent energy resolution, see 
sect. 5'4. 
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kilograms to several hundreds, depending on the proposal. 

Unfortunately, the (3/3 isotopes are not always abundant in nature, requiring enrich- 
ment in order to obtain large, concentrated masses. It was argued in [92] that ^^^Xe 
would be a particularly favorable isotope to use, since it permits target masses of 1 ton 
or more and low-background experimental techniques. Isotope enrichment appears rela- 
tively easier from a technical point of view (and therefore, cheaper) for ^"^^Xe. Experimen- 
tal proposals involving both liquid scintillator detectors (see KamLAND-Zen, sect. 6'6), 
as well as liquid-phase or gaseous-phase TPCs (sec EXO and NEXT, sects. 6'3 and 6'7, 
respectively), are using or planning to use ^'^^Xe. 

Liquid scintillator proposals permit in principle to reach large f3/30v isotope masses 
with other isotopes as well. Given its very favorable phase space, the SNO-I- Collaboration 
(see sect. 6'8) will dissolve a neodymium salt in the liquid scintillator. Scalability to large 
i50]\j(j masses will ultimately depend on the feasibility to enrich neodymium, most likely 
a difficult enterprise. 

High-resolution calorimeters, that is germanium diodes and bolometers, are compact 
detectors and might be therefore, in principle, scalable to large masses. In germanium 
experiments such as GERDA (see sect. 6'4), /3/3 isotope masses of the order of 100 kg 
appear feasible, but it might be difficult, in particular for economical reasons, to go much 
beyond that. The GUORE bolometers (see sect. 6'2) plan to use ^'^°Te. This isotope 
has the highest natural isotopic abundance among the commonly-considered [3f3 emitters 
(34%). The need for isotope enrichment is therefore less important in this case, and 
/3f3 masses of the order of hundreds of kilograms appear within reach of new-generation 
experiments. 

5'4. Energy resolution. - Together with a large isotope mass, good energy resolution is 
a necessary (but not sufficient!) requirement for the ultimate (3[3Qu experiment. It is the 
only protection against the intrinsic I3j32v background, and improves the signal-to-noise 
ratio in the region of interest around Qpfj. 

The detectors for (3(3 searches that have achieved the best energy resolution so far 
are the germanium diodes and the bolometers. In germanium detectors the energy is 
measured via ionization (creation of electron- hole pairs in the semiconductor) . An energy 
resolution as low as 0.1% FWHM at Qpp has been obtained [96]. Partly thanks to their 
superior energy resolution, germanium diodes have dominated the /3/30i^ searches so far 
(see sect. 6T). The GERDA and MAJORANA proposals are based on germanium diodes, 
see sect. 6'4. In bolometers the energy is measured by detecting a temperature rise in 
crystals with very small specific heat. Several bolometric crystals have been proposed and 
tested for /3/3 searches, with tellurite (Te02) being the favored one due to its reasonable 
mechanical and thermal properties, and the natural high content (28% in mass) of the 
^'^'^Te /3/3 isotope. An energy resolution of about 0.2% FWHM at Qpp has been reached 
using Te02 crystals [66]. The GUORE proposal uses this technique (see sect. 6'2). 

5'5. Backgrounds. - Double beta decay experiments are mostly about suppressing 
backgrounds. As we have seen already, the mere presence of background in the region 
of interest around Qpp changes the regime of the mpp sensitivity from a {Mpp ■ t)~^/^ 
dependence to [Mfjp ■ t)~'^/^. 

The natural radioactivity of detector components is often the main background in 
[3j3Qv experiments. Even though the half-lives of the natural decay chains are comparable 
to the age of the Universe, they are very short compared to the half-life sensitivity of the 
new-generation [3j3Qi' experiments. Therefore, even traces of these nuclides can become a 
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Fig. 19. - Decay chains of uranium (left) and thorium (right) [97] . 



significant background. The decays of ^os^pj ^^^^ ■^^'^Bi are particularly pernicious, given 
the high Q- values of these reactions, therefore polluting the energy region of interest of 
most /3/3 emitters. These isotopes are produced as by-products of the natural thorium 
and uranium decay chains (see fig. 19), and they are present at some level in all materials. 
Careful selection of material and purification is mandatory for all PfiQv experiments. The 
new-generation experiments are being fabricated from amazingly radiopure components, 
some with activities as low as 1 /^Bq/kg or less. 

Radon gas, either ^^^Rn or ^^°Rn, is also a worry for most experiments. These 
isotopes, present in the natural decay chains, diffuse easily through many materials, 
infiltrating the detectors sensitive region. Their daughters tend to be charged and stick to 
surfaces. Many experiments eliminate radon from the detector surroundings by flushing 
pure nitrogen gas. Also, some laboratories have installed radon-traps in the air circulation 
system. 

In addition to internal backgrounds coming from radioactive impurities in detector 
components themselves, there are external backgrounds originated outside the detec- 
tor. Those backgrounds can be in principle suppressed by placing the detector at an 
underground location and by enclosing it into a shielding system. 

In general, the depth requirement for a PI3Qv experiment varies according to the de- 
tector technology. A very efficient shielding and additional detection signatures such 
as topological information can compensate the benefits of a very deep location. Figure 
20 shows the depth (and corresponding cosmic ray muon flux) of several underground 
facilities currently available to host physics experiments around the world. The deep- 
est laboratory is SNOLAB (Canada), an expansion of the existing facilities constructed 
for the Sudbury Neutrino Observatory (SNO). The SNO-f- experiment (sect. 6'8), and 
perhaps also EXO (sect. 6'3), will be located there. Other deep laboratories include 
SUSEL (USA) and LSM (France), which wiU host the demonstrators for the MA JO- 
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Fig. 20. - Total muon flux as a function of the equivalent vertical depth for a flat overburden. The 
empirical parametrization, shown as a dashed line, is taken from [99]. The fluxes measured at 
various underground sites currently available to host physics experiments are taken from [99, 98]. 
Facilities shown in brown, blue and red are located in Europe, America and Asia, respectively. 
The full names of the facilities shown in the figure are: Waste Isolation Pilot Plant (WIPP), 
Laboratorio Subterraneo de Canfranc (LSC), Soudan Underground Laboratory (SUL), Kamioka 
Observatory (Kamioka), Boulby Palmer Laboratory (Boulby), Laboratorio Nazionale del Gran 
Sasso (LNGS), Laboratoire Souterrain de Modane (LSM), Sanford Underground Science and 
Engineering Laboratory (SUSEL), Baksan Neutrino Observatory (BNO), SNOLAB. 



RANA (sect. 6'5) and SuperNEMO (sect. 6'9) experiments, respectively. In addition to 
depth, other important factors characterizing the underground sites include the size of 
the excavated halls and the services provided to the experiments. The size is an impor- 
tant factor to take into account especially for experimental proposals at the ton-scale, 
given that some of them {e.g., the SuperNEMO experiment) need large volumes. For 
a recent review of the currently available underground facilities around the world, see 
reference [98]. 

At the depths of underground laboratories, muons (and neutrinos) are the only sur- 
viving radiation from cosmic rays. However, their interactions can produce high-energy 
secondaries such as neutrons or electromagnetic showers. 

Charged backgrounds (such as muons) can be easily eliminated using a veto system. 
Neutrons, on the other hand are often a more serious problem. They can have sizable 
penetrating power, impinging on the detector materials and activating them through 
large AA transitions in nuclei, ultimately resulting in radioactive nuclides. Cosmogenic 
activation is, of course, more severe on surface. Therefore, for experiments using materials 
that can get activated (like germanium-based experiments), underground fabrication and 
storage of the detector components is essential. The detectors can be shielded against 



THE SEARCH FOR NEUTRINOLESS DOUBLE BETA DECAY 



41 



neutrons with layers of hydrogenous materiaL 

Radioactive decays in the rock of the underground cavern resuh in a 7-ray flux that 
can interact in the detector producing background. Dense (high Z), radiopure materials 
such as lead and copper are used as shielding to suppress this background. Water, being 
inexpensive and easy to purify, is also a good alternative for shielding against 7-rays. 

Finally, very massive detectors such as liquid-scintillator calorimeters suffer from an 
irreducible external background: the solar neutrino flux. 

In the design of a shielding system against external backgrounds, a graded shielding 
principle is followed: the thickness of a shield component does not need to reduce the flux 
below the contribution of the next inner component, with the innermost shield component 
selected to be the radiopurest. 

The lowest background rate (expressed in terms of background events per unit en- 
ergy, /3/3 isotope mass and exposure time) in a (U/SOv experiment so far was achieved by a 
tracker-calo experiment, NEMO-3 (see sect. 6'1). In this experimental approach, foils of 
the Pf3 source are surrounded by a tracking detector that provides a direct detection of 
the two electron tracks emitted in the decay. The topological reconstruction of the events 
provides a powerful active handle to reject backgrounds, together with relatively radiop- 
ure detectors (see fig. 21). The NEMO-3 experiment measured a background rate of a few 
times 10""^ counts/ (keV • kg^^ • year) [100]. Time projection chambers using xenon in gas 
phase, as proposed by the NEXT experiment (see sect. 6'7), provide also some topological 
information that can be used to reject backgrounds. High-resolution calorimeters, such 
as germanium diodes and bolometers, have so far achieved somewhat worse background 
rates, in the neighborhood of 10~^ counts/(keV • kg^^ • year) [1, 101, 66]. The goal 
of the new-generation experiments is typically to reach 10""^ counts/(keV • kg^^ • year) 
background levels, and sometimes significantly better than that. This may require, in 
some cases, activities for detector components of 1 /iBq/kg or less! A recent example 
witnessing the challenges to build and characterize such extremely radiopure detectors 
is provided by the EXO Collaboration. In this case, a thorough and systematic study of 
trace radioactive impurities in a large variety of parts and materials required to construct 
the EXO-200 detector has been carried out and documented in detail [102]. 

Another handle to suppress non-/3/32;y backgrounds is daughter ion tagging. This has 
been proposed, and is actively being pursued, for the EXO xenon-based detector. In this 
case, the f3 pOiy decay is ^"^^Xe — >■ ^^^Ba+"'' -I- 2e~. The -'^•^^Ba+'*" ion rapidly captures an 
electron, resulting in ^"^^Ba"'' which is stable in xenon. The ^^^Ba+ ions can be identified 
via atomic spectroscopy, exciting them with a blue laser and observing the resulting 
red light (sec fig. 22). Daughter ion tagging is undoubtedly very challenging from the 
technical point of view, but the payoff would be huge if the R&D were to be successful. 

5'6. Detection efficiency. - Neutrinoless double beta decay events are extremely rare, 
if present at all, thus a high detection efficiency is an important requirement for a (3(3 
experiment. Equation (52) clearly indicates that the detector design should prioritize 
a high detection efficiency. To obtain the same increase in mpp sensitivity obtained by 
doubling the efficiency, the mass would have to be increased by a factor of 4, assuming 
the same background. 

In general, the simpler the detection scheme, the higher the detection efficiency. For 
instance, pure calorimetric approaches such as germanium diodes or bolometers have 
detection efficiencies in excess of 80%. This is to be contrasted with experiments per- 
forming, for example, particle tracking, which will typically cause a significant efficiency 
loss. Also, homogeneous detectors, where the source material is the detection medium. 
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Fig. 21. - Top (left) and side (right) view of a reconstructed event selected from NEMO-3 
data with a two electron energy sum of 2812 keV [34]. 
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Fig. 22. - Simplified energy level diagram of the barium ion, and the wavelength in vacuum for 
the transitions. 



provide in principle higher efficiency than the separate-source approach. This is due 
to a number of reasons, including geometric acceptance, absorption in the /3/3 source, 
backscattering of electrons, and the tracking requirement. On the other hand, some 
relatively dense homogeneous detectors use some of the /3/3 mass close to the detector 
borders effectively for self-shielding, paying it with some efficiency loss. 

6. A selection of new-generation experimental proposals 

During decades the search for double beta decay was a rather marginal activity carried 
out with geochemical techniques. It was not until 1987 that the /3/32i' was directly ob- 
served in the laboratory. In the 1990's, the field was dominated by germanium detectors, 
devices characterized by superb energy resolution and a high efficiency. 

After the positive results of neutrino oscillation experiments, the field has gone through 
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a revolution. The community is preparing a rich and varied new generation of exper- 
iments that should explore ultimately the inverted-hierarchy region of neutrino masses 
(see fig. 10). This will require a multi-ton experiment. It seems prudent to build, as a 
first step, an experiment containing about 100 kg of isotope that can be expanded at a 
later time. However, the scalability will not be possible for all experimental techniques. 

In this section, after summarizing the results of past experiments, some proposals for 
the new generation are described. This discussion does not pretend to be exhaustive, 
but focused on the pros and cons of the different techniques. Finally, the sensitivity to 
of the proposals is evaluated. 

6'1. Past experiments. - For almost half a century the only evidence of the existence of 
double beta decay came from geochemical methods consisting in measuring the concen- 
trations of the stable daughter isotopes {Z + 2, A), produced over geologic times 10^ 
years). An excess of the daughter isotope over its natural concentration is interpreted 
as evidence for /3/3 decay (either /3f32v or fS/SQi/, since the method cannot distinguish 
between them). 

The first direct measurement of /3/32i^, in ^^Sc, did not happen until 1987 [22]. It 
was done using a fairly large 1 m'^) time projection chamber, the well-known Irvine 
TPC. The source, 14 g of 97% enriched ^^Se, was deposited on a thin Mylar foil forming 
the central electrode of the chamber. The trajectories of the electrons emitted from the 
source foil were recorded by the TPC and analyzed to infer their energy and kinematic 
characteristics. Since this initial detection, the two-neutrino mode has been directly 
observed for 8 isotopes in several experiments (see table II and ref. [23] for further 
details). 

The most restricting limits to date in the search for P/SOv were obtained with ger- 
manium detectors. The Heidelberg-Moscow (HM) experiment [1] searched for the /3/30j^ 
decay of ""^Ge using five high-purity Ge semiconductor detectors enriched to 86% in ''^Ge. 
The experiment ran in the Laboratori Nazionali del Gran Sasso (LNGS), Italy, from 1990 
to 2003, totaling an exposure of 71.7 kg-year. The background rate reached by the ex- 
periment in the region was (0.19 ± 0.01) counts/(keV • kg • year), or, in units of 
(3/3 emitter mass, 0.22 ± 0.01 counts/(keV • kgpfj ■ year). Pulse shape discrimination 
(PSD) was used in a subset of the data (35.5 kg-year) to separate single-site events, like 
PPOv decays, from multi-site events, like 7 interactions, resulting in a background rate 
of (0.06 ± 0.01) counts/(keV • kg • year), or (0.07 ± 0.01) counts/(keV • \gpp ■ year). A 
lower limit on the /3/30j/ half-life oiT°J^{^^Ge) > 1.9 x lO^^ years (90% CL) was obtained 

[!]• 

A subset of the collaboration re- analyzed the data claiming evidence for ''^Ge P/SOv 
decay [68]. The latest publication by this group reports a 6a evidence for /3f3Qv and a 
half-life measurement of T^^^ = (2.23^0 3^) x 10^^ years [70], corresponding to m/^p — 

(0.30^0 03) eV according to the central value of the PMR nuclear matrix clement for 
^^Ge given in sect. 4'3. This claim sparked an intense debate in the community, and at 
the moment no consensus exists about its validity (sec, for example, ref. [69]). 

The International Germanium Experiment (IGEX) [71] also searched for fS/SQi/ using 
enriched germanium crystals. It ran in the Homestake gold mine (USA), the Canfranc 
Underground Laboratory (Spain) and the Baksan Neutrino Observatory (Russia) from 
1991 to 2000, accumulating a total exposure of 8.87 kg-year. It reached a sensitivity 
similar to that of Heidelberg-Moscow, but not enough to disprove the claim. The lowest 
background rate reached by the IGEX experiment was 0.26 (0.10) counts/ (keV • kg • year) 
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without (with) pulse shape discrimination for a 8.87 (4.65) kg-year total exposure [101], 
corresponding to 0.30 (0.12) counts/(kcV • kg^^ • year) per unit /3/3 emitter mass. 

The Cuoricino experiment, an array of 62 Te02 bolometric crystals, ran for five years 
in Gran Sasso searching for (^/30iy in -'^'^^Te. It reached a sensitivity to to^/j comparable to 
that of the HM experiment, but it cannot disprove the claim due to the uncertainties in 
the nuclear matrix elements. The average background rate for the 5x5x5 cm"^ Cuoricino 
crystals, computed in a 60 keV wide region centered around Qfj/s, was 0.161 ± 0.006 
counts/(keV • kg • year) [103], corresponding to 0.58 ± 0.02 counts/(keV • kg^^ • year) 
per unit /3/3 emitter mass. The average FWHM energy resolution in all crystals was 
6.3 ± 2.5 keV at 2615 keV [103]. 

The lowest levels of background so far were achieved by the NEM03 experiment [100]: 
a few times 10"'^ counts/(keV • kg • year). This detector represents the state of the art 
of separate-source /3/3 experiments. Reconstruction of the electron tracks emerging from 
the source provided a powerful signature to discriminate signal from background. The 
NEM03 experiment ran from 2003 to 2010 at the Modane Underground Laboratory 
(LSM), in France. The detector, of cylindrical shape, had 20 segments of thin source 
planes, with a total area of 20 m^, supporting about 10 kg of source material. The sources 
were within a drift chamber, for tracking, surrounded by plastic scintillator blocks, for 
calorimctry. A solenoid generated a magnetic field of 25 Gauss which allowed the mea- 
surement of the tracks electric charge sign. The detector was shielded against external 
gammas by 18 cm of low-background iron. Fast neutrons from the laboratory environ- 
ment were suppressed by an external shield of water, and by wood and polyethylene 
plates. The air in the experimental area was constantly flushed, and processed through a 
radon-free purification system embedding the detector volume. In addition to searching 
for /3/30t^, NEM03 very successfully served as a "/3/32i/ factory", providing precise /3/32z/ 
half-life measurements for seven /3/3 isotopes, see table II. Apart from representing the 
ultimate background for /3/30i^ searches , an accurate measurement of j3l32u in several 
nuclides is also important as input to NME calculations. 

6'2. CUORE. - The Cryogenic Underground Observatory for Rare Events (CUORE) 
[104] has been designed following the successful experience of the MiDBD [44] and Cuori- 
cino [66] ^'^'^Te experiments, where for the first time arrays of bolometers were used to 
search for j3(3 decay. 

CUORE will be placed in the hall A of the Gran Sasso Underground Laboratory and 
will consist of a system of 988 bolometers, each being a crystal of Te02 of 5 x 5 x 5 
cm"^, arranged in 19 vertical towers consisting of 13 layers of 4 crystals each. The four 
crystals are held between two copper frames joined by copper columns. PTFE pieces are 
inserted between the copper and Te02, as a heat impedance and to clamp the crystals. 
There is a few mm gap between crystals with no material between them. A system 
of lead shields will be hosted inside the cryostat close to the detectors, to shield them 
from environmental radioactivity and from radioactive contaminations originating in the 
dilution unit located above the detector and in the dewar structure [105]. A 6 cm thick 
roman lead shield surrounds the detector array on its sides, while a 30 cm thick layer of 
low-activity lead is placed above the detector. The ^^°Pb activity of the roman lead was 
measured to be less than 4 mBq/kg. A sketch of the detector is shown in fig. 23. 

The total mass of the detector will be 741 kg for a ^'^°Te mass of 206 kg. The energy 
released in a single particle interaction within the crystal is measurable as a change 
in temperature by Neutron Transmutation Doped (NTD) germanium thermistors. The 
measured energy resolution is ^ 5 keV FWHM at the /3/3 {Qv) transition energy 2.53 




Fig. 23. - The CUORE detector and cryostat. In yellow, the 19 towers of bolometers; the 
lavender volumes are the lead shielding. 



McV). 

The CUORE bolometers will operate at temperatures between 10 and 15 mK. A 
challenging '^He/'^He dilution refrigerator, with a cooling power of 3 mW at 120 mK, has 
been designed on purpose and is under construction. 

A single tower of CUORE, CUORE-0, is presently under construction and will begin 
operations within 2011. It will be hosted in the old Cuoricino dilution refrigerator, placed 
in the hall A of the LNGS. CUORE-0 is a real test of the CUORE assembly chain and 
procedure, will directly test the level of backgrounds of the CUORE setup and improve 
the Cuoricino sensitivity on jSpQv. 

The CUORE setup will possibly allow in the long term for powerful upgrades. An 
obvious, though expensive, possibility is to substitute the natural tellurium bolometers 
with enriched ^^°Te units (provided that the enrichment procedure can keep the internal 
backgrounds very low). A more sophisticated option is to use scintillating crystals con- 
taining interesting double beta emitters [106]. The contemporary read-out of scintillation 
light and thermal signal could indeed allow for a dramatic reduction of the background 
rate and a better characterization of the signal. An array of ZnSe scintillating bolome- 
ters, LUCIFER, has been recently proposed as a prototype experiment exploring the 
performances of such an approach [107] (see sect. 6'10). 

6-3. EXO. - The Enriched Xenon Observatory [108] wiU search for /3^0i/ in ^^exe. The 
ultimate goal of the Collaboration is the development of the barium tagging for a multi- 
ton xenon-based detector, which would lead to a virtually background-free experiment. 
Prior to that, the Collaboration has built the EXO-200 detector, a '-^200-kg liquid xenon 
(enriched to 80.6% in ""^Xe) time projection chamber that detects both scintillation and 
ionization. 

The fiducial volume of the chamber, 44 cm in length, is divided in two halves by a 
central cathode (see fig. 24, left). Ionization charges created in the xenon by charged 
particles drift under the influence of an electric flcld towards the two ends of the cham- 
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Fig. 24. - Left: one half of the EXO chamber, viewed from the cathode plane. Right; the 
chamber attached to the cryostat door, as viewed from the bottom of the APD plane. The legs 
contain the readout cabling and are also the conduits for xenon circulation. 



ber. There, the charge is collected by a pair of crossed wire planes which measure its 
amplitude and transverse coordinates. Each end of the chamber includes also an array 
of avalanche photodiodes (APDs) to detect the 178-nm scintillation light produced by 
primary interactions. The sides of the chamber are covered with teflon sheets that act as 
VUV reflectors, improving the light collection. The simultaneous measurement of both 
the ionization charge and scintillation light of the event may in principle aUow to reach a 
detector energy resolution as low as 3.3% FWHM at the ^^^^Xe Q- value, for a sufficiently 
intense drift electric field [109]. 

The xenon is held inside a thin copper vessel immersed in a cryofiuid that also shields 
the experiment from external radioactive backgrounds. The HFE heat-transfer fiuid is 
stored in a vacuum-insulated low-activity copper cryostat. The cryostat is surrounded 
on all sides by 25 cm of low-activity lead. The entire assembly is surrounded by a radon- 
free tent and housed in a class 100 clean room, the exterior of which is instrumented on 
five sides with plastic scintillator panels for vetoing cosmic rays with 95.9% efficiency. 
The detector is located 2150 feet underground for an overburden of 1585 meters water 
equivalent, at WIPP (Waste Isolation Pilot Plant), in the United States. 

The EXO-200 TPC was installed in its cryostat in early 2010. The detector was 
first filled with natural (unenriched) xenon in late 2010. The data collected during 
these engineering runs were used to make a first assessment of the performance of the 
detector, and to perform a first round of calibrations. Low-background running with 
enriched xenon started in the spring of 2011. The EXO Collaboration announced the 
observation of the /?/32i/ mode of ^'^^Xe in August 2011. Prior to this measurement, 
reported in table II, the f3[i2v had been observed in all other important PfiQv candidate 
nuclei except ^'^^Xe. The measured /3/32zy half-life is significantly lower than previously 
reported lower limits [67, 110]. 

To identify the daughter barium, several methods are under study, including single- 
ion fluorescence, resonant ionization spectroscopy (RIS), and mass spectroscopy. Single 
ion fluorescence is a highly sensitive and highly selective method to observe a barium ion 
while held under vacuum in a RF trap. In this technique, the Ba"*" ion is rapidly cycled 
from its 6^Si/2 ground state to its 6^Pi/2 excited state by illuminating it with lasers of 
the appropriate wavelength (493 nm and 650 nm), see fig. 22. As the electronic state 




Fig. 25. - Sketch of the GERDA experiment. The germanium arrays can be seen inside the 
copper cryostat, and this one placed inside the cylindrical water tank. 



changes, the laser photons are scattered in all directions, and the scattered light can be 
easily detected by a photo-multiplier tube. EXO has achieved good single barium ion 
identification with this technique, even in the presence of low pressure xenon and helium 
gas mixtures. However, this technique also requires that the barium ion be retrieved 
from the TPC volume, transported to the RF trap, released, and trapped, while not 
altering its chemical or ionization state. Resonant ionization spectroscopy, on the other 
hand, is a technique which allows single barium ions to be observed without requiring 
a vacuum ion trap. In RIS, barium ions are desorbed from the surface of a transport 
probe, and subsequently resonantly ionized under illumination by 554 nm and 390 nm 
lasers. The ionized barium can then be observed with a Channeltron electron multiplier. 
Initial tests with the RIS technique have successfully identified barium being desorbed 
from the probe tip, so this technique is promising. Other avenues of research include 
barium identification within xenon ice, and barium extraction from a high pressure gas 
TPC using gas nozzles. 

6'4. GERDA. - The GERmanium Detector Array (GERDA) experiment [111], located 
in Hall A of the Laboratori Nazionali del Gran Sasso (LNGS), will make use of naked 
Ge detectors immersed in a large cryostat of ultra-pure LAr. 

The Ge detectors are organized in strings (2-5 detectors) and mounted in special 
low-mass (^ 80 g) holders made of ultra-pure copper and PTFE. The array of strings 
is contained in a vacuum insulated stainless steel cryostat of 4.2 m diameter and 8.9 m 
height. A copper shield covers the inner cylindrical shell of the cryostat with a maximum 
thickness of 6 cm. The cryostat is placed in a water tank, of 10 m diameter and 9.4 m 
height, serving as a gamma and neutron shield. It will be also used as a veto against cos- 
mic rays thanks to its instrumentation with 66 photomultipliers, with good efficiency in 
detecting the Cherenkov light. The cosmic muon veto is reinforced by plastic scintillator 
panels on top of the detector, for a surface of about 20 m^. A drawing of the detector 
and shielding is shown in fig. 25. 

In its first phase, GERDA-1, eight fully refurbished germanium diodes (17.7 kg total 
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active mass, 86% isotopic enrichment in ^^Ge) from the previous Heidelberg-Moscow and 
IGEX experiments will be used. In the subsequent step, GERDA-2, new diodes will be 
used for a total active mass of 35.4 kg. These new diodes will be p-type Broad Energy 
(BEGe) detectors [112, 96], allowing for a better discrimination of backgrounds thanks 
to a sophisticated pulse shape discrimination. 

The experiment started commissioning runs in June 2010 using natural Ge, low- 
background, detectors, refurbished from the Genius- TF experiment. These commission- 
ing runs focused on the investigation of the background sources (most notably on the un- 
expectedly large contribution from ^^Ar-^'^K in the LAr) and on ''^Ar background mitiga- 
tion strategics, see sect. 6' 12. Data taking for the GERDA-1 physics run will start in the 
next months. The background level of the natural Ge setup was measured to be 0.06±0.02 
counts/ (keV • kg • year) (corresponding to 0.07±0.02 counts/(keV • kg^^j • year)), consis- 
tent with early indications from the first string of enriched Gc detectors deployed [113]. 
This rate, obtained without using pulse shape information, is a factor of 3-4 lower than 
the HM and IGEX measured ones (see sect. 6T), but still about a factor of 6 higher 
than the GERDA-1 goal. The reason for this higher than expected background rate is at 
present not fully understood. The goals of GERDA-2 are to start data taking in about 
two years, with about twice the isotope mass of GERDA-1, and with a background level 
of 0.001 counts/(keV • kg • year) (or 0.0012 counts/(keV • kg^^ • year)). 

In the very long term a third phase of the experiment, GERDA-3, is foreseen to 
make use of about 1 ton of ''^Ge target material together with a further reduction of 
background. Such an effort, common with the MAJORANA project (see below), would 
be feasible only in a word-wide collaboration, and provided that the GERDA approach 
could demonstrate to be the best candidate technology to push the double beta decay 
sensitivity below the inverted hierarchy mass threshold (about 30 meV). 

6'5. MAJORANA. - The MAJORANA Collaboration is following a more classic ap- 
proach than GERDA in the design of a germanium-based experiment [114]. The Ge 
detectors will be mounted in a string-like arrangement in ultra-pure vacuum cryostat 
made from radiopure copper. The cryostat will be surrounded by a passive shielding of 
Cu and Pb, and an active muon veto. 

The Collaboration is building a demonstrator module, to be placed at the Deep Un- 
derground Science and Engineering Laboratory (DUSEL) in the United States, with 
about 20 kg of natural BEGe detectors. The goal is to demonstrate a background rate 
of about 4 counts per tonne and per year in the 4-keV wide region of interest [115]. The 
demonstrator is expected to operate with enriched detectors in 2013. 

6'6. KamLAND-Zen. - The KamLAND-Zen experiment [116] will search for /3/30z^ in 
^^^Xe using enriched xenon dissolved in liquid scintillator. This will allow a calorimetric 
measurement of the [3(3 electrons, as first proposed in [117]. Xenon is relatively easy to 
dissolve (with a mass fraction of more than 3% being possible) and also easy to extract 
from the scintillator. 

The major modification to the existing KamLAND detector [118] was the construction 
of an inner, very radiopure (of order 3x 10""'^^ g/g of ^^^U and '^'^^Th) and very transparent 
balloon to hold the dissolved xenon. This balloon, 1.58 m in radius, is placed at the center 
of the KamLAND active volume as shown in fig. 26. 

The KamLAND-Zen experiment plans to dissolve 389 kg of -'^'^^Xe in the liquid scin- 
tillator of KamLAND in the first phase of the experiment, and up to 1 ton in a projected 
second phase. 
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Fig. 26. - Sketch of the KamLAND-Zen detector. The ball-on containing the dissolved xenon 
(purple) hangs in the center of the active volume. 



The proven resolution (from the previous operation of the KamLAND experiment) is 
16% FWHM at 1 MeV. The main sources of expected background are the f3l32i> tail, ^^^Bi 
impurities in the scintillator or in the balloon, ^"^C generated in the scintillator by cosmic 
rays, and ®B solar neutrinos. The expected background rate in the region of interest is 2 x 
10~^ counts/(keV • kg • year) [119], corresponding to 2.2x10"^ counts/(keV • kg^g^ • year). 

At the time of writing this report, the mini-balloon installation into the KamLAND 
detector has been completed, and detector commissioning is ongoing. Physics data-taking 
with the xenon- loaded liquid scintillator is expected to start in the fall of 2011. 

67. NEXT. - The Neutrino Experiment with a Xenon TPC (NEXT) [120] will search 
for PPOv in ^^^Xe using a 100-kg high-pressure gaseous xenon (HPXe) time projection 
chamber. Such a detector can provide both good energy resolution and event topological 
information for background rejection [121]. 

Double beta decay events leave a distinctive topological signature in HPXe: a ion- 
ization track, of about 30 cm long at 10 bar, tortuous due to multiple scattering, and 
with larger energy depositions at both ends (see fig. 27). The Gotthard experiment [122], 
consisting in a small xenon TPC (5.3 kg of xenon, 68% enrichment in ^^^Xe) operated 
at 5 bar, proved the effectiveness of such a signature to reject background, achieving a 
background rate of only ^ 0.01 counts/(keV • kg • year). 

The design of NEXT is optimized for energy resolution (better than 1% FWHM 
at Qpfj) by using proportional electroluminescent (EL) amplification of the ionization 
signal. The detection process is as follows. Particles interacting in the HPXe transfer 
their energy to the medium through ionization and excitation. The excitation energy is 
manifested in the prompt emission of VUV (~178 nm) scintillation light. The ionization 
tracks (positive ions and free electrons) left behind by the particle are prevented from 
recombination by a strong electric field (0.5-1.0 kV/cm). Negative charge carriers drift 
toward the TPC anode, entering a region, defined by two highly-transparent meshes. 
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Fig. 27. - Simulation of a PPOu track in gaseous xenon at 10 bar [120]. 



with an even more intense electric field (3.5 kV/cm/bar). There, further VUV photons 
are generated isotropically by electroluminescence. Therefore, both scintillation and 
ionization produce an optical signal, to be detected with a sparse plane of PMTs located 
behind the cathode. The detection of the primary scintillation light constitutes the 
start-of-event (to), whereas the detection of EL light provides an energy measurement. 
Electroluminescent light provides tracking as well, since it is detected also a few mm away 
from production at the anode plane, via a dense array (1 cm pitch) of 1-mm^ SiPMs. 

The NEXT detector will operate at 10 bar, with xenon enriched at 90% in the ^^^Xe 
isotope. At that pressure the 100 kg mass of xenon results in a volume of ~2.5 m'^. 

The major benefits of the NEXT 100 proposal are its high background rejection factor, 
resulting in an expected background rate of 2 x 10~* counts/(keV • kg^/j • year), and the 
fact that xenon is relatively easy (cheap) to enrich and obtain in large quantities. 

The NEXT Collaboration expects to commission the detector at the end of 2013. The 
experiment plans to start its physics run in the second half of 2014. 

6-8. SNO+. - SNO+ [123] is the follow-up of the successful SNO experiment [124], 
located at SNOLAB, in Canada. It re- uses the existing equipment of the detector 
(acrylic vessel, photomultipliers and their support structure, electronics and the light 
water shield) replacing the heavy water by ^780 tonnes of liquid scintillator (linear 
alkylbenzene, LAB). 

The physics program of the SNO+ detector includes measurements of low energy solar 
neutrinos and /3(30i/ searches using ^^°Nd. In order to do that, the liquid scintillator will 
be loaded with a neodymium salt, resulting in about 50 kg of ^^°Nd. This isotope has 
the second highest endpoint, 3.37 MeV, and the fastest predicted neutrinoless double 
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Fig. 28. - One of the candidate configurations for the SNO+ acrylic vessel anchor system. The 
acrylic vessel is shown in grey, and the anchor system in yellow. The outer sphere made of 
triangles is the PMT support structure. 



beta decay rate due to its large phase space factor, see fig. 12. The high endpoint is 
above most radioactive backgrounds, such as radon, and this is a significant advantage. 
However, enrichment of this isotope seems difficult. 

The energy resolution of the SNO+ detector is estimated to be 6.5% FWHM at 3.4 
MeV. External backgrounds can be rejected with a relatively tight fiducial volume selec- 
tion, cutting however about 50% of the signal. The most important sources of background 
are expected to be ^os^pj impurities in the scintillator, the irreducible background from 

solar neutrinos and P[i2v events from ^^°Nd. Assuming the radiopurity levels for 
the liquid scintiUator achieved by BOREXINO lO-^"^ g/g of ^o^Tl) [125], simulations 
predict a background rate of ~ 10^^ counts/(keV • ^gpp ■ year) [126]. 

The SNO-I- experiment is expected to start commissioning in the spring of 2013 with 
pure liquid scintillator, to be followed by the Nd-loaded liquid scintillator phase. Given 
that the LAB liquid scintillator is about 15% less dense than the surrounding light water, 
one of the major technical challenges of the SNO-I- upgrade is the design of a hold-down 
system for the acrylic vessel using a net of radiopure ropes, see fig. 28. 

6'9. SuperNEMO . - This proposed new installment of the NEMO detectors series 
consists of up to 20 tracker-calo modules, each one containing a thin foil of about 5 kg 
of /3/3-decaying material, probably ^^Se, although other isotopes such as ^^"Nd or ^^Ca 
are also under consideration. 

A sketch of a SuperNEMO module can be seen in fig. 29. The source foil, 3 meters high 
and 4.5 meters long, with a surface density of about 40 mg/cm^, is placed in the center 
of a tracking chamber with overall dimensions of 4 m height, 5 m length and 1 m width. 
Nine planes of drift cells operating in Geiger mode and a magnetic field of 25 Gauss allow 
to reconstruct the trajectory and charge of particles crossing the chamber. A calorimeter 



52 J.J. GOMEZ-CADENAS, J. MARTIN-ALBO, M. MEZZETTO, F. MONRABAL and M. SOREL 




Fig. 29. ~ A SuperNEMO module. The source foil (not shown) is placed in the center of a 
tracking volume consisting of drift cells operating in Geiger mode. The tracking volume is 
surrounded by calorimetry consisting of scintillator blocks connected to PMTs (grey). The 
support frame is shown in red. 



consisting of blocks of plastic scintillator coupled to low-activity PMTs surrounds the 
tracking chamber on four sides. Its granularity allovifs the energy of individual particles 
to be measured. 

The physics case of SuperNEMO relics on several significant improvements over the 
NEMO-3 detector performance [127]. The energy resolution is expected to be 7% FWHM 
at 1 MeV, a factor of 2 better than in NEMO-3. Such a resolution has been attained 
with a 28 cm hexagonal PVT scintillator directly coupled to a 8-inch PMT [128]. The 
detection efficiency of SuperNEMO is estimated by means of simulation to be about 30%, 
almost a factor of 2 better than in NEMO-3. As far as the backgrounds are concerned, 
SuperNEMO goals require an impressive improvement in the purification (both chemical 
and via distillation methods) of the source foils. In particular, ^^'*Bi and ^os^j contami- 
nation in ^^Se foils are to be reduced by factors of 50 and 170, respectively. A dedicated 
setup, the BiPo detector, installed in the Laboratorio Subterraneo de Canfranc (LSC), 
will measure the radiopurity of the foils in order to make sure that the required levels 
are achieved. Finally, in order to decrease radon gas levels in the tracking chamber down 
to negligible levels (<0.15 mBq/m'^ ), a reduction of at least a factor of 40 with respect 
to NEMO-3 is needed. 

The first SuperNEMO module, called the demonstrator, will be the first step from 
R&D to construction with the aims to demonstrate the feasibility of large scale mass pro- 
duction, to measure the backgrounds (especially from radon emanation), and to finalize 
the detector design. The demonstrator will be installed in the space previously occupied 
by the NEMO-3 detector at the Modane Underground Laboratory. 

The current plans of the SuperNEMO Collaboration for the following: (a) demon- 
strator construction, 2010-2012; (b) demonstrator physics run start-up, 2013; and (c) 
full detector construction start-up, 2014. 



6'10. Other proposals. - 
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CANDLES. This project [129] proposes the use of CaF2 scintillating crystals to search 
for f3(3Qv in ^®Ca. The crystals would be immersed in liquid scintillator providing shield- 
ing and an active veto against external backgrounds. Among the /3/3 isotopes, ^^Ca has 
the highest Q-value, 4.27 MeV. This places the signal well above the energy region of 
the natural radioactive processes. Unfortunately, the natural abundance of the isotope 
is only 0.187% and enrichment seems complicated. Therefore, many tons of crystals are 
needed for a competitive new-generation experiment. 

COBRA. The COBRA experiment [130, 131] is exploring the potentials of Cadmium 
Zinc Telluride (CdZnTe) room-temperature semiconductor detectors for PfiQv searches. 
Out of the several /3/3 candidate isotopes in CdZnTe, COBRA is focusing on ^'^°Te, 
because of its natural abundance, and ^^®Cd, because of its high Q-value of 2.8 MeV. 
Activities are split in two main directions: (a) the identification of the main background 
components in a setup of 64 commercial 1-cm'^ CdZnTe diodes located at LNGS; and 
(b) the development of pixelized devices that would allow to reduce the background by 
particle identification. 

DCBA. The Drift Chamber Beta-ray Analyzer [132] is a magnetized tracker (drift 
chambers) that can reconstruct the trajectories of charged particles emitted from a /3/3 
source foil. The momentum and kinetic energy are derived from the track curvature in the 
magnetic field. A prototype, DCBA-T2, has shown energy resolution of about 150 keV 
(FWHM) at 1 MeV, and the main source of background (^-'^^Bi) has been identified. A 
new apparatus, DCBA-T3, with a more intense magnetic field is now under construction 
at KEK. 

LUCIFER. The idea of LUCIFER [133, 107] is to join the bolometric technique pro- 
posed for the CUORE experiment with the bolometric light detection technique used in 
cryogenic dark matter experiments. Preliminary tests on several jS^Qv detectors have 
clearly demonstrated the background rejection capabilities that arise from the simulta- 
neous, independent, double readout (heat and scintillation). LUCIFER will consist of 
an array of ZnSe crystals operated at 20 mK. The proof of principle with about 10 kg of 
enriched Se is foreseen for 2014. 

MOON. The MOON detector [134] is a stack of multi-layer modules, each one con- 
sisting of a scintillator plate for measuring energy and time, two thin detector layers 
for position and particle identification, and a thin /3/3 source film interleaved between 
them. At present, Nal(Tl) scintillators are considered as the candidates for the scintil- 
lator plates. Energy resolution around 3% FWHM at 3 MeV has been achieved during 
the R&D phase. For position-sensitive detectors, possible candidates are multi-wire pro- 
portional chambers (MWPCs) and Si-strip detectors. 

XMASS. XMASS [135, 136] is a multi-purpose liquid xenon scintillator. Although 
optimized for dark matter searches, it will also investigate neutrinoless double beta decay 
and solar neutrinos. The detector, with about 800 kg of xenon, was installed in the 
Kamioka mine (Japan) in the fall of 2010. The excellent self-shielding capabilities of the 
liquid xenon will be used to define a virtually background-free inner volume. 
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6'11. Sensitivity of new- generation experiments. - In this section we try to assess 
the physics case of the new-generation double beta experiments described above(^^). 
We quote the experimental sensitivities to to^^, assuming the standard light Majorana 
neutrino exchange as the dominant f3/30v mechanism. To perform this risky exercise we 
make use of the physics-motivated ranges for the NME values described in sect. 4'3, and 
of the set of experimental parameters summarized in table IV. A discussion motivating 
our choice of parameters in table IV is given in sect. 6T2. 



(^^) For the sensitivity computation, we restrict ourselves to experiments that involve at least a 
few kg of (3/3 emitter mass, that are approved, and that have been granted a significant financial 
support. 



Table IV. - Basic parameters for the different double beta experiments: I3j3 emitter mass Mpp, j3pQv efficiency e, FWHM energy resolution 
AE, and background rate c per unit energy, /3/3 isotope mass and time. The last column indicates the number of background events within the 
ROI, and is the product of the Mpp, /S.E and c columns. Comparison of different approaches is very difficult directly from the numbers in the 
table, but this information is fundamental to compute their sensitivity. 



Experiment 




e 


AE 


c 


Bgr/ROI 








(keV) 


(10~^ counts/(keV • kg/j^ ■ year)) 


(cts/yr) 


EXO-200 


141 


0.34 


100 


0.78-5 


11-71 


GERDA-1 


15.2 


0.95 


4.2 


12-70 


0.77-4.5 


GERDA-2 


30.4 


0.84 


2 


1.2-7 


0.07-0.43 


CUORE-0 


10.9 


0.83 


5 


180-390 


9.8-21.3 


CUORE 


206 


0.83 


5 


36-130 


37.1-134 


KamLAND-Zen 


357 


0.61 


250 


0.22-1.8 


19.6-161 


MAJORANA Demonstrator 


17.2 


0.85 


2 


1.2-12 


0.04-0.41 


SNO+ 


44 


0.50 


220 


9-70 


87-680 


NEXT 


89.2 


0.33 


18 


0.2-1 


0.32-1.6 


SuperNEMO Demonstrator 


7 


0.28 


130 


0.6-6 


0.55-5.5 
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Although different experimental aspects are relevant from the point of view of the 
feasibility of an experiment, the sensitivity can be computed using only a few parameters 
— namely, effective mass of the isotope and background rate in the energy Region Of 
Interest (ROI) — that can be extracted from the fundamental numbers in the design of 
the experiments. 

The sensitivity is calculated based on the Feldman-Cousins method [137] for construct- 
ing confidence intervals, following the prescription given in [92]. For each experiment, 
we define a ROI centered at the /3/3 decay Q-value and extending for one FWHM of 
energy resolution, and we compute the experimental sensitivity at 90% confidence level. 
We take into account the effect of the FWHM selection (corresponding to 76% efficiency 
assuming gaussian resolution) as a multiplicative factor to the experimental efficiency 
reported in table IV. 

Our prescription assumes a counting experiment in the ROI with a perfectly known 
background rate. In other words, in our sensitivity computation, we neglect systematic 
uncertainties(^^) and any energy shape information that may be present in the energy 
distribution of events. Systematic uncertainties may possibly affect the parameters listed 
above, especially the knowledge of the backgrounds, and deteriorate the sensitivity. On 
the other hand, use of additional information beyond the overall count rate of /3/30i/ 
candidates within the ROI may yield some sensitivity improvement. While important, 
both effects would be extremely difficult to incorporate in such a sensitivity comparison, 
given that most new-generation experiments discussed here have not even started their 
commissioning phase, yet. 

In fig. 30 we show the m/^p sensitivities at 90% CL for the new-generation proposals 
discussed above, assuming a 5 years exposure for all of them. The colored rectangles 
reflect the uncertainty coming from the PMR nuclear matrix elements, see fig. 16. For 
each experiment, two rectangles are shown, corresponding to the most optimistic and 
most pessimistic background rate expectations given in tab. IV, respectively. For most 
experiments, the two rectangles overlap. For each experiment, the thin solid line at the 
lowest mpp values is meant to give an idea of what can be gained by increased exposures. 
The line represents the m^^ sensitivity for the most optimistic NME values (within the 
PMR range of fig. 16) and for the most optimistic background rate expectations (within 
the ranges of tab. IV) for a 10 years exposure. For reference, we also show the mass 
range of the standing KK claim of evidence for /3/30iy in ^^Ge [70] (with a central value 
of 300 mcV) and the mass region as predicted under the inverted hierarchy hypothesis 
(between 17 and 52 meV, see fig. 10). The sensitivity of the various proposals after a 5 
years exposure, both in terms oiT^J^ and mpj^, is also shown in tab. V, where the central 
value of the PMR interval for the nuclear matrix elements and the optimistic background 
rate expectations in tab. IV have been assumed. 

The first thing to remark from fig. 30 is that the KK claim should be unambiguously 
solved by several new-generation proposals using different isotopes. If our assumptions 
are correct, this will certainly be the case for '^^Ge (GERDA, MAJORANA), "Oxe 
(CUORE) and "6Xc (EXO-200, KaniLAND-Zcn, NEXT), and possibly also for ^^^Se 
(SuperNEMO) and i^°Nd (SNO-f-). Muhi-isotope determination of /S/SOi' may therefore 
become a real possibility within this decade. In this respect, it is important to note that 
GERDA and MAJORANA are the only experiments among those in fig. 30 using the 



(^^) With one exception for the SNO+ and KamLAND-Zen experiments, see sect. 6'12, where 
an attempt has been made to account for systematic effects affecting energy reconstruction. 
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EXO200 GERDA1 CUORE-0 KaniLAND-Zen SNO+ GERDA2 CUORE NEXT SuperNemo D. Majorana D. 

Experiments 

Fig. 30. " Sensitivity of the different experiments to the neutrino effective mass m^/j computed 
assuming a 5 years exposure, the PMR intervals for the nuclear matrix elements (see sect. 4'3), 
and for both the optimistic and pessimistic experimental parameters of table IV. A statistical 
90% CL is computed according to the Feldman-Cousins method [137] , assuming a signal region 
of one FWHM and the corresponding efficiency. For each experiment, the sensitivities for the 
two experimental parameter sets are drawn as overlapping rectangles. A sensitivity line corre- 
sponding to a 10 years exposure, and to the most optimistic NME and experimental parameter 
set, is also shown. The upper grey region represents the KK claim [70] while the lower grey 
region represents the inverted hierarchy region (see fig. 10). 



same isotope of the HM experiment (^•^), and should therefore be able to provide the only 
truly model-independent confirmations of this claim. 

On the other hand, several experiments appear to have a very good chance to reach 
a sensitivity of 100 meV or better, in particular CUORE, KamLAND-Zen, NEXT and 
EXO-200. In our most optimistic scenario concerning NME values and experimental 
performances, this target may also be reached by GERDA during its second phase. 
Given our uncertainties, we cannot predict vifhich, among these 4-5 different experimental 
proposals, will provide the best to^^ sensitivity after a 5 years exposure. To this end, a 
better knowledge of the actual (as opposed to expected) values for the background rates, 
of the systematic uncertainties affecting the measurement, and of the NME values would 
be necessary for all experiments. 

From fig. 30, our expectation is that it will be almost impossible for the new- generation 
experiments discussed here to discover P^Qv after 5 years of exposure if the neutrino mass 



(^^) Not only: we have seen that GERDA in its first phase re-uses the same detectors as in the 
HM experiment. 
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Table V. - Sensitivity of the experiments at 90% CL after a 5 years exposure, both in terms of 
half-life Ty2 and in terms of neutrino effective mass mpp. These values are obtained from the 
optimistic background rate assumptions in tab. IV. The conversion from TIJ2 mpp assumes 
the central value of the PMR interval for the nuclear matrix elements. 



Experiment 


T°j2 (years) 


mpii (meV) 


CUORE-0 


8.67 X 10^^ 


203 


CUORE 


8.86 X lO^'^ 


63 


GERDA-1 


4.49 X 10^^ 


252 


GERDA-2 


1.37 X 10^*^ 


121 


EXO200 


8.20 X IQ^^ 


82 


NEXT 


9.13 X lO^'^ 


78 


KamLAND-Zen 


1.32 X lO^*^ 


65 


SNO+ 


5.38 X 10^'' 


182 


SuperNEMO Demonstrator 


9.15 X lO^'^ 


258 


MAJORANA Demonstrator 


7.19 X lO^'' 


258 



spectrum is hierarchical (miight — 0) rather than degenerate, since essentially no over- 
lap exists between the experimental sensitivities and the 17-52 meV inverted hierarchy 
region. Not only larger exposures, but also new (better) experimental proposals, would 
be needed to fully probe this mass region. As mentioned above, we expect experiments 
using i36xe (EXO-200, KamLAND-Zen, NEXT) to provide for the first time during this 
decade a comparable or better sensitivity than ^"^Ge (GERDA, MAJORANA) and ^^''Te 
(CUORE) experiments, which dominated the field over the past two decades. In per- 
spective, if the low-background expectations of new-generation ^^^Xe proposals will be 
confirmed during this decade, ^'^^Xe would be a particularly favorable isotope to use also 
in the longer term. This is because scalability to isotope masses in the 1-10 ton range 
are in this case more realistic than with any other isotope. 

Figure 30 represents our main result for the physics case comparison of new-generation 
PpOv experiments. In this figure, a single, physics-motivated, NME uncertainty band is 
used for each isotope, following our discussion in sect. 4'3. For completeness, we have 
also repeated the same exercise for several other NME values or ranges, one for each the- 
oretical framework considered in [85]. The result is shown in fig. 31. As mentioned above, 
the spread of the corresponding predictions most likely overestimates the theoreti- 
cal uncertainty in the rnpp sensitivity conversion. For figure clarity, only nuclear 
theory assumptions are varied in fig. 31, while the detector performance parameters are 
fixed to the most optimistic values of tab. IV. As can be seen in fig. 30, our assumed un- 
certainties in the detector performance parameters would yield m^p sensitivity changes 
of about the same size as the nuclear theory variations shown in fig. 31. 

6T2. Validity of sensitivity assumptions . - The reliability of the sensitivity estimates 
given above critically depend on how realistic our choice of detector performance indi- 
cators is. In this section, we discuss how we have chosen the parameters reported in 
tab. IV. This discussion is mostly intended for the expert reader wishing to indepen- 
dently assess our choices, and to use his/her own judgment to modify them accordingly. 
Given that the largest uncertainty affecting the ultimate mpp sensitivity of a proposal is 
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EXO- 200 GERDA-1 CUORE- KamLAND- Zen SNO+ GERDA-2 CUORE NEXT SuperNEMO D. Majorana D. 

Fig. 31. " Sensitivity of the different experiments to the neutrino effective mass mpp computed 
assuming the optimistic experimental parameters of table IV. A statistical 90% CL is computed 
according to the Feldman-Cousins method [137], assuming a signal region of one FWHM and 
the corresponding efficiency. Five different frameworks for NME calculations are considered, 
following reference [85], and are drawn as overlapping rectangles. The upper grey region rep- 
resents the KK claim [70] while the lower grey region represents the inverted hierarchy region 
(see fig. 10). 



almost always related the achievable background rates, we decide to quote a background 
rate range. For all other indicators, a single number rather than a range is used, see 
tab. IV 

The EXO-200 TPC is filled with about 175 kg of liquid xenon enriched to 80.6% in 
the isotope ^'^^Xe [24], corresponding to a /3/3 mass of about 141 kg^^. For the /3/30i^ 
efficiency around Q/3/3, we take the efficiency assumed by the EXO Collaboration for their 
/3(32i/ analysis [24], corresponding to e = 0.34 above the 720 keV analysis threshold. The 
inefficiencies are dominated by the fiducial volume cut, keeping 63 out of 175 kg of liquid 
xenon [24] (0.36 efficiency). A 6.3% inefficiency introduced by vetoing /3-a coincidences 
[24] has also been considered in our efficiency assumption. In [24], the collaboration 
measured an energy resolution of as/E = 4.5% at 2615 keV for the EXO-200 detector. 
This value was obtained using a 376 V/cm drift field and ionization signals only. An 
improvement of up to a factor of 2.5 could be achieved with higher (1-4 kV/cm) drift 
fields and combining ionization with scintillation information, see [109]. As a result, 
we assume a nearly-nominal, 4.2% FWHM energy resolution at Qpi3, corresponding to 
about 100 kcV. For our background rate lower limit, we consider the collaboration's 
goal of 20 radioactive background events per year in a ±2a interval around the Qj3i3 
endpoint, for a a/E = 1.6% energy resolution at 2.5 MeV, a detector mass of 200 kg 
and a 80% enrichment in ^"^^Xc [138]. From these numbers, we obtain a background 
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rate of c = 0.78 x 10"'^ counts/(keV • kg^^ • year). This nominal background rate 
prediction still remains to be updated based on real EXO-200 data. As a worst-case 
background rate scenario, we take the rate that has already been achieved: 4 x 10"'^ 
counts/ (keV • kg • year), see [139]. This background level was obtained without full 
lead shielding, radon exclusion tent, radon trap or full 3-dim reconstruction, and might 
therefore be improved in the future [139]. This number corresponds to c = 5 x 10^'^ 
counts/(keV • kg^^ • year) per unit /3/3 mass. 

GERDA-1 will use eight refurbished ''^Ge diodes from the Heidelberg-Moscow and 
IGEX experiments, for a total active mass of 17.66 kg and 0.86 isotopic enrichment in 
™Ge [140], corresponding to a /3/3 mass of 15.2 kg/3/3. As for the /3/30z^ efficiency, the 
actual value will ultimately depend on analysis details that are unknown at the moment, 
for example whether the collaboration will rely on pulse shape discrimination already 
in phase I to further reduce multi-site energy deposition events. In the absence of an 
updated number, we assume e = 0.95 as originally quoted by the collaboration in [141]. 
The FWHM energy resolution for GERDA-1 diodes was measured to be between 3.6 
and 6.0 keV at the 2615 keV gamma ray line from ^os^j^ [113]. Taking the central value 
of this interval (4.8 keV) and extrapolating to the ^^Ge value (2.039 MeV), we 
estimate AE = 4.2 keV. The optimistic background rate scenario is assumed to be 
the collaboration's goal of 0.01 counts/(keV • kg • year) [141], corresponding to 0.012 
counts/(keV • kg^^ • year) per unit /3/3 mass. GERDA started commissioning in mid- 
2010, and data obtained since then can be used to estimate a worst-case scenario for 
the achievable background rates. The best background rate measured with a string of 3 
natural germanium detectors, refurbished from the Genius-TF experiment, is 0.06 ±0.02 
counts/(keV • kg • year) [113]. Since mid-2011, the first enriched germanium detectors 
have been deployed on a second string arm, using the best detector configuration found 
so far. Preliminary data from the enriched germanium detectors indicate a background 
rate that is compatible with the one found with the natural germanium diodes [113]. As 
a consequence, we assume 0.06 counts/(keV • kg • year) as upper limit for the GERDA-1 
background rate, translating into c = 0.07 counts/(keV • kg^^ • year) per unit P(3 mass. 

For the phase II of the experiment, the GERDA Collaboration purchased 37.5 kg 
of germanium with an isotopic abundance in ^^Ge of 0.86. The material has already 
been purified into 35.4 kg of 6N germanium [142], corresponding to a /3/3 mass of 30.4 
kg^^. In order to reduce backgrounds, both sophisticated pulse-shape discrimination 
(PSD) techniques and additional instrumentation for the LAr veto are likely to be used 
by GERDA-2. We assume an overall /3/?0i^ efficiency of e = 0.84, given by the product 
of the 0.86 efficiency for a PSD cut reported in [143], times the 0.973 efficiency for a 
LAr veto cut reported in [144]. Compared to phase I detectors, a significantly improved 
FWHM energy resolution of 2 keV a.t Qpp has been measured for the BEGe detectors 
to be used by GERDA-2, see [96]. The lower limit on the background rate is taken to 
be the collaboration's goal of 0.001 counts/(keV • kg • year) [141], corresponding to c = 
0.0012 counts/(keV • kg^^ • year). Again, we use prehminary results from the GERDA 
commissioning runs to estimate an upper limit on the background rate. Commissioning 
data indicate that /? decays of ^^K, that is in turn produced positively charged by the 
^^Ar decay within the LAr veto, can contribute very significantly to the Qpp background 
rate. This ^^K background can be most easily quantified by measuring the 1525 keV 
gamma ray line. On the one hand, the collaboration estimated via simulations that a 
background rate at Q^p of up to 1.7 x 10"'^ counts/ (keV • kg • year) can be obtained by a 
homogeneous distribution of ^^K around the detectors, for a 43.9 /iBq/kg contamination 
in '*^K [145]. On the other hand, a 1525 kcV line about 20 times more intense than 
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expected was observed during the first commissioning run. A significant fraction of this 
enhancement has been understood as due to a inhomogcneous ^^K distribution caused by 
the field fines drifting the positively-charged ^^K ions to the detector surface. To prevent 
the ^^K ions to reach the detector surfaces, the collaboration deployed a copper shield 
called the "mini-shroud". The additional shroud reduced the counts at the 1525 kcV line 
and at Qp/s by a factor of 4-5, and a preliminary measurement of the *^Ar specific activity 
of about 160 ^Bq/kg was obtained for almost field- free runs [140]. This measured value, 
combined with the simulation result of 1.7 x 10^"^ counts/(keV • kg • year) for a 43.9 
/iBq/kg contamination in ^^K, motivates our upper limit background rate assumption of 
about 0.006 counts/(keV • kg • year), or c = 0.007 counts/(kcV • kg/j^ • year) per unit 
/3(3 mass. 

CUORE-0 wiU make use of a single CUORE-like tower with 39 kg of natural Te02 
crystals and with an isotopic abundance of 0.34167 of ^'^'^Te [103], corresponding to a (3/3 
mass of 10.9 kg^^. The overall signal efficiency has been estimated to be about 0.83, as 
obtained for the big crystals of the CUORICINO experiment [66] . The main inefficiency 
source is the "physical" inefficiency due to beta particles escaping the detector and ra- 
diative processes. The expected FWHM energy resolution of the CUORE detectors is 
AE ~ 5 keV at the (3/301' transition energy [103]. As lower limit on the background rate, 
we assume 0.05 counts/(keV • kg • year) from the 2.615 keV gamma ray multi-Compton 
events coming from the irreducible ^"^^Th contamination of the CUORICINO cryostat, 
as done in [103]. This rate corresponds to c = 0.18 counts/(keV • kg^^ • year) per unit 
(3(3 mass. As upper limit on the background rate, again we follow [103] and assume 0.11 
counts/ (keV • kg • year), translating into c = 0.39 counts/ (keV • kg^^ • year) per unit (3(3 
mass. This latter number includes an additional background contribution from scaling 
the CUORICINO background in the conservative case of a factor of 2 improvement in 
^■^^U and ^"^^Th contamination of the copper and crystal surfaces. Such surface con- 
tamination results in degraded alphas that may mimic the (3(3 signal. This factor of 2 
improvement is largely motivated by the background rates measured in the Three Towers 
Test (TTT), allowing to estimate the surface contamination of the copper detector hold- 
ers, responsible for a large (50 ± 20% [103]) fraction of the CUORICINO backgrounds 
[146]. For such tests, crystals were dismounted from the CUORICINO detector, and 
repolished on the surfaces. Three different types of copper cleaning were tested. The 
copper treatment procedure selected by the collaboration proved to be able to reduce 
the copper surface contamination by at least a factor of 2 compared to CUORICINO 



The CUORE total active mass will be 741 kg for the 988 envisaged cubic detectors 
[103]. As for CUORE-0, the natural Te02 crystals will have an isotopic abundance of 
0.34167 of ^30Te [103], resulting in about 206 kg^^j of "°Te. The overafi signal efficiency 
and the FWHM energy resolution at Qpp are assumed to be equal to the CUORE-0 
figures given above, e ~ 0.83 and AE ~ 5 keV, respectively. The main improvements 
in background rate reduction compared to CUORE-0 are assumed to come from a much 
cleaner cryostat and from the different detector geometry. Under these assumptions, the 
CUORE background will be dominated by copper and crystal surface contaminations. 
As lower limit on the background rate, the 0.01 counts/(keV • kg • year) goal of the col- 
laboration is assumed [103, 147], corresponding to c = 0.036 counts/(keV • kg^^ • year). 
As upper limit, an extrapolation of the currently available measurements [146] on cop- 
per and crystal surface contaminations to the CUORE geometry results in a 0.035 
counts/(keV • kg • year) upper limit [148], resulting into c = 0.13 counts/ (keV • kg^^ • year) 
per unit f3(3 mass. 



[103, 146]. 
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During the phase I of the experiment, KamLAND-Zen will use 389 kg of xenon [119], 
enriched to 0.917 [149] in the isotope ^■^^Xc, for a ^■^^Xe total mass of 357 kg^/j. Concern- 
ing f3(30i' efficiency, the collaboration will make use of a fiducial volume cut to suppress 
backgrounds that are external to the mini-balloon, for example due to ^■^^U/^'^^Th con- 
taminations in the outer liquid scintillator or the environmental gamma ray background. 
Depending on the ^'^^U/^'^^Th contamination levels of the mini-balloon materials them- 
selves (Nylon film, supporting ropes and pipes), the collaboration is also considering a 
tighter fiducial volume cut to suppress backgrounds originating from the mini-balloon. In 
the following, we assume that a fiducial volume definition placed 3 cr,- inward with respect 
to the mini-balloon surface will be used, where (Tr is the radial position reconstruction 
accuracy. Assuming a,. = 12.5 cm/ y^E{MeV) [119], we estimate a corresponding fidu- 
cial volume efficiency of 0.61 for the 158 cm radius mini-balloon deployed. No other 
sources of inefficiency are considered in our estimate. Concerning energy reconstruc- 
tion, the same performance measured in the previous phase of experiment is expected 
for KamLAND-Zen, given that the xenon-loaded liquid scintillator within the balloon 
has the same optical properties (light yield, transparency) as the KamLAND scintilla- 
tor outside the balloon. An energy resolution of as/E = 6.8%/ ^E(McV) is assumed 
[119], scaling to 250 keV FWHM energy resolution at Qpp- For the background rate 
lower limit, we take the latest collaboration's expectations, obtained via simulations. 
In [119], 19.5 background events per year are expected. The three largest background 
sources are expected to be I3(32i/ events entering the ROI, ^^^Bi events from the mini- 
balloon materials, and (to a lesser extent) "'^^C events produced through cosmic muon 
spallation. Compared to previous estimates [149], a shorter (3/321/ half-life as measured 
in [24] is assumed {T1/2 ~ 2 x 10^^ years), together with a less radiopure mini-balloon 
(^'^^U/^'^^Th concentrations of 3 x 10~^^ g/g)- This estimate results in a background rate 
of c = 0.22 X 10~'^ counts/(keV • kg^^ • year). Factors that may potentially result in a 
highcr-than-cxpected background rate are a non-perfect knowledge of the reconstructed 
energy spectrum of (3/321/ events spilling over the ROI, a higher background contribution 
from mini-balloon materials, and a worse tagging of ^^^Bi and ^"C backgrounds (esti- 
mated tagging efficiencies of 66% and 90%, respectively). It is, however, rather difficult 
to quantitatively estimate what a "pessimistic" background rate might be observed, at 
this stage. As in the SNO-f- discussion below, we assume (to a large extent in a arbi- 
trary fashion) a 8 times higlier-than-expected background rate as upper limit. We note 
that more information to revisit the KamLAND-Zen background model should become 
available soon, given that KamLAND-Zen data-taking has started. 

The MAJORANA demonstrator will contain 40 kg of germanium, of which at least 20 
kg and up to 30 kg will be enriched to 86% in ^^Ge [115]. We follow the collaboration's 
baseline and assume 20 kg of enriched germanium [150], corresponding to 17.2 kg^^ 
of "^^Ge. A total (3(301/ efficiency of 71% is estimated in [151], accounting for detector 
granularity, pulse-shape analysis (PSA), single-site time correlation (SSTC) and energy 
cuts. A 5% loss due to edge effects and lost bremsstrahlung is also included in this 
number. We divide by the quoted 84% efficiency of the energy cut alone to estimate a 85% 
efficiency in tab. IV. For energy reconstruction, the 2 keV FWHM resolution measured 
at Qpp in BEGe detectors [96] is assumed, as for GERDA-2. For the background rate, we 
assume as lower limit the number quoted by the collaboration: 0.004 counts per year and 
kg in a 4 kcV-widc ROI [115], corresponding to c = 1.2 x 10~^ counts/(kcV • kg^^ • year). 
Main backgrounds responsible for this rate are expected to be prompt cosmogenics, and 
208rpj^2i4gj contaminants in the cryostat and in the copper/lead shield [152]. A non- 
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negligible contribution from ^®Ge contaminants in the enriched germanium crystals is also 
expected. The material purity specifications are extremely challenging, with contaminant 
goals at the < 0.1(< 0.3) /iBq/kg level in ^os^pj (^^''^Bi) for the electroformed copper to 
be used for detector mounts, cryostat and the inner copper shield, and about one order 
of magnitude worse for the commercial high-purity copper and lead to be used for the 
outer shield [151, 152]. Purities within one order of magnitude of such goals have already 
been demonstrated, see [151]. Wc therefore assume a factor of 10 worse background rate 
than nominal as upper limit: c ~ 12 x 10^^ counts/(keV • kg/j^ • year). 

The SNO-I- experiment will make use of 780 tonnes of liquid scintillator [153], with 
natural neodymium loading at the 0.1% (w/w) [126]. Given the 5.6% isotopic abundance 
of ^^°Nd, this concentration corresponds to 44 kg/j^ of /3/3 emitter mass. The experiment 
is expected to have a fiducial volume cut rejecting about 50% of the active volume [126]. 
No other sources of inefficiency are considered in our estimate. For energy reconstruction, 
we assume a light output of 400 NHit/MeV [126], or 1347 NHit at Q^^, where NHit is 
the number of PMT hits. This number is highly sensitive to the amount of neodymium 
concentration, and corresponds to 0.1% (w/w) loading. Such a light output implies a 
FWHM resolution of about 220 keV at Qpf;. For the background rate, of the order of 100 
background events per kton of liquid scintillator and per year are expected via simulations 
in a 200 keV energy window around Qpp [126]. This background level translates into 
a rate of c = 0.009 counts/(keV • kg^^ • year), which we take as lower limit. The 
three dominant backgrounds are expected to be ®B solar neutrinos, ^os^p^ decays and 
/3/32i/ events. The background level due to ^B solar neutrinos is well-known. The ^ospp 
background assumes a level of ^'^^Th impurities at the level measured by the BOREXINO 
experiment, 8.3 x lO"""^® g/g. The collaboration expects energy reconstruction systematic 
effects, such as non-gaussian resolution tails, on the shape of the /3/32j/ energy spectrum 
near the ROI to affect the sensitivity the most, see for example the study in [126]. 
Preliminary studies conducted by the collaboration indicate that a m^p sensitivity within 
a factor of 5/3 worse than the purely statistical sensitivity can be preserved including 
this systematic effect [154]. In the large background approximation, which is valid for the 
SNO-I- experiment, this systematic effect would therefore be equivalent to a background 
rate increase of up to a factor of (5/3)^ ~ 8, resulting in a background rate upper limit 
of c = 0.07 counts/(keV • kg^^ • year). 

The NEXT detector will contain 99.14 kg of pressurized xenon gas in the TPC fidu- 
cial volume region, enriched to 0.90 in the isotope ^■^^Xe [155], corresponding to 89.2 
kg^^ in mass. A /3/30i/ efficiency of 0.25 has been estimated via simulations in [155], 
accounting for inefficiencies in reconstructing the /3/3 track, and in imposing energy and 
topology cuts to suppress backgrounds. Given that the inefficiency introduced by the 
energy within ROI requirement is separately accounted for in our analysis, we assume an 
efficiency of e = 0.25/0.76 = 0.33 in tab. IV. Preliminary energy reconstruction measure- 
ments obtained with a kg-scale prototype yield a 4.6% FWHM resolution at the 59.4 keV 
full energy peak for a ^^^Am calibration source [155]. This energy resolution extrapolates 
to 0.72% FWHM resolution at (5/3/3, or about 18 keV. As lower limit on the background 
rate, we take the collaboration's estimate of 0.2 x 10^'^ counts/(keV • kg/3^ • year), dom- 
inated by the 238u^232rpj^ contamination of the titanium pressure vessel, conservatively 
assumed to be at the level of 200 /j,Bq/kg for each isotope, see [155]. Such radiopurity 
assumptions arc based upon the measured upper limits for the clean titanium used in 
the LUX cryostat [156]. As upper limit for the background rate, we take an independent 
(and more pessimistic) measurement of titanium radiopurity, at the level of 1.2 ± 0.4 
(0.6 ± 0.3) mBq/kg for ^ssy (232^^-]^^ using the Gator low-background counting facility 
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at LNGS [157, 158]. Contaminants in these amounts would translate into a background 
rate of about 10"'^ counts/(kcV • kg^/j • year) at Q/j/j. Also, the collaboration has started 
its own radiopurity R&D campaign at LSC, with the goal of refining these assumptions 
in the near future. 

For SuperNEMO, we assume a 7 kg mass in the isotope ^^Se, to be installed in the 
demonstrator module [127]. A /3/30j^ efficiency of 0.28 has been estimated for a Su- 
perNEMO module in a detailed study [159], accounting for acceptance, reconstruction 
efficiency and event selection efficiency. We assume this number in our estimates, which 
is in fact quite similar to the collaboration's goal of e = 0.30 [128]. Calorimeter R&D 
efforts achieved a 7.7%/^E(McV) FWHM energy resolution using a PVT scintillator 
directly coupled to a 8 inch, high QE, Hamamatsu PMT, see [128]. This energy reso- 
lution measurement extrapolates to 4.4% (or about 130 keV) FWHM at Qpp. For the 
background rate, we take c = 6 x 10"'^ counts/ (keV • kg^/j • year) as worst-case scenario. 
This number comes from the preliminary result on the fi^Qv search in ^^Se with the 
NEMO-3 detector [160]: 14 events were observed (in agreement with background expec- 
tations) in the [2.6-3.2] MeV energy ROI, after 4.5 years of data-taking and 0.93 kg/3/3 
of source foil. The backgrounds are dominated by ^^^Bi/^°®T1 contamination of the foils 
(measured to be 530 ± 180 ^Bq/kg and 340 ± 50 /xBq/kg in ®^Se, respectively, see [100]) 
and radon concentration in the tracking volume (6.46 ± 0.02 mBq/m'^ for the phase-2 
of the experiment, see [100]). As background rate lower limit, we assume a factor of 10 
improvement in both ^'^^Tl/^^'*Bi radiopurity of the foils and in radon concentration in 
the tracker: c = 0.6 x 10~^ counts/(kcV • kg/3/3 • year). 

7. Conclusions 

Double beta decay is a rare nuclear transition in which a nucleus with Z protons decays 
into a nucleus with Z+2 protons and the same mass number A. The standard double beta 
decay mode, the two-neutrino mode /3/32j/, is a second-order weak transition producing 
two electrons and two antineutrinos, and is allowed in the Standard Model. Despite 
being a very slow process, it has been observed in a variety of even-even nuclei. The 
neutrinoless double beta decay mode /SpOv is instead a hypothetical process producing 
two electrons and no neutrinos. While several other processes have been investigated, 
[ij3Qv is the most promising probe we have in hand to test lepton number violation. Also, 
a positive result in the search for PpQv would unavoidably indicate that neutrinos are 
Majorana particles, that is truly neutral particles. 

The experimental exploration of PPQv has a long history, see fig. 32. The first search 
was performed in 1948, where the half-life for jSfiQv in ^^^Sn was constrained to be 
longer than 3 x 10^^ years [161]. Half a century later, in 2001, the Heidelberg-Moscow 
Collaboration reported a half- life limit about ten orders of magnitude more stringent, 
using ^^Ge as /3/3 emitter: T^j2 > -'-•^ ^ ^^^^ years [1]. It is, also, a history plagued by 
frequent claimed discoveries (see, for example, [162]) that have been later disproved by 
subsequent experiments. This observation alone reflects how difficult it is to search for 

There is at present a diverse and healthy competition among a variety of experimental 
techniques to establish themselves as the best approach for fifiQv searches. The (ifiQv 
field is now witnessing a golden age in terms of experimental efforts. Why is that? Some 
reasons have been present all along during the era of P^Qv exploration: 

• We have a fairly good idea of what to look for. While several mechanisms have 
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Fig. 32. - Seventy years of direct PfiOv searches in perspective. Existing limits (sliown in blue) 
are taken from [161]. The sensitivity of new-generation proposals (shown in red) is based upon 
this review, see sect. 6. 



been proposed to drive /3/30;/, in most of them the two decay electrons arc the only 
light particles emitted, therefore carrying most of the available energy. This can 
be contrasted with proton decay searches, where it is less clear which decay mode 
should be the focus of experimental investigation. 

• It is common belief that there is still ample room for improvement with respect 
to the most sensitive /3(30i/ searches performed to date, as can be guessed by the 
trend in fig. 32. 

There are, however, additional reasons that are applicable to the present era: 

• Probably the most important reason has to do with the discovery of neutrino os- 
cillations over the past two decades, implying that neutrinos are massive particles. 
If one assumes, as it is customarily done, that light Majorana neutrino exchange 
is the dominant contribution to /J/JOi^, there is a direct link between a measurable 
PPOv rate on the one hand, and the absolute scale of neutrino masses scale and neu- 
trino oscillations phenomenology on the other. In this context, one can also study 
what is the actual value of neutrino masses, and whether the neutrino mass spec- 
trum exhibits some particular features (such as a hierarchical or a quasi-degenerate 
structure), via PfSQv. 



• Searching for jSpOv is well motivated on theoretical grounds. On the one hand, there 
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is no fundamental reason why total lepton number should be conserved. On the 
other hand, Majorana neutrinos provide natural explanations for both the smallness 
of neutrino masses and the baryon asymmetry of the Universe. As a consequence, 
theoretical prejudice in favor of Majorana neutrinos has gained widespread consen- 
sus. 

• As it is well known, there is a 6 ct evidence for /3/30z^ in ^^Ge reported by part of 
the Heidelberg-Moscow Collaboration, T^"^^ = (2.23l!^;3i) x lO^^ years [70]. It is 
also well known that this claim is highly controversial [69]. Consensus exists that 
the issue can only be definitely settled by new, and more sensitive, experiments. 

The mapping of observed /3/30;/ rates into neutrino mass constraints not only requires 
assuming the standard /3/30i/ interpretation in terms of light Majorana neutrino exchange. 
It also requires precise nuclear physics knowledge, which can be factorized into the so- 
called nuclear matrix elements (NMEs). These NMEs cannot be measured, and need 
to be separately calculated for each /?/3 emitting isotope under consideration. Several 
calculations exist. While they share common ingredients, calculations differ in their 
treatment of nuclear structure. We argue that about a 20-30% NME uncertainty exists 
for converting rates into neutrino masses. 

In this review, the different experimental aspects affecting P(30iy searches were ex- 
tensively discussed. The requirements are often conflicting, and no new-generation ex- 
perimental proposal is capable of optimizing all of them. Should we concentrate on 
approaches offering huge event rates, such as KamLAND-Zen? Are superior energy reso- 
lution techniques, such as GERDA or CUORE, the best way to go? Should background 
suppression focus more on radiopurity control, as in the EXO or KamLAND-Zen cases, 
or on powerful signal-background discrimination techniques, as in the SuperNEMO or 
NEXT approaches? We made an attempt at a quantitative comparison of the physics 
case of selected new-generation experimental approaches, which is summarized in fig. 30. 

What about the longer-term future? How far can we go in fHpOiy exploration? Nobody 
knows for sure. What we do know is that such future (3f30i/ searches will unavoidably need 
to involve experiments at the ton or multi-ton scale in /3/3 isotope mass. The diversity 
of experimental approaches wc arc currently witnessing will not be viable at that scale, 
and 2 or 3 approaches (most likely based on different isotopes) are going to be retained 
at most. However, an extrapolation of the trends from the past and the present may 
offer some qualitative clues. Figure 32 seems to point to "asymptotic" limits of I3(30i' 
half-life explorations in the 10^^ years range. If the light Majorana neutrino exchange 
mechanism is realized in Nature, this would correspond to effective Majorana masses at 
the few meV scale. As can be appreciated in fig. 10, such ultimate /3/30z^ sensitivities 
would give us good chances to detect fSpOiy regardless of the value of the neutrino mass 
and mixing parameters. 

An unambiguous detection, either in the current-generation efforts starting now or in 
the longer-term future, would open up an even more exciting era for /3/30i/ searches, with 
the objective to actually understand what is the physics mechanism that is responsible 
for this elusive process. 
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